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ABSTRACT(U) 

This  document  represents  a  user  manual  for  the  Turbo  Pascal  /GEM  software  interface  that 
has  been  developled  recently  within  Guided  Weapons  Division.  The  interface  has  been 
developed  primarily  to  enable  Turbo  Pascal  programmer  to  replicate  their  softcopy 
graphical  output  in  high  quality  hardcopy  form.  The  output  from  the  interface  is  a  GEM 
file  that  can  be  edited  as  required  within  either  GEM  Draw  or  GEM  Artline,  prior  to 
generation  of  a  PostScript  file.  A  number  of  mathematical  procedures  have  been  included 
in  order  to  extend  the  capabilities  of  Turbo  Pascal,  together  with  a  useful  library  of  curve 
fitting  procedures.  The  interface  also  incorporates  a  form  filling  technique  for  versatile  data 
entry  atany  point  within  a  Turbo  Pascal  program.  A  detailed  description  of  all  procedures 
and  functions  is  included  together  with  the  code  and  graphical  output  from  a  number  of  test 
programs  designed  to  illustrate  the  capabilities  of  the  interface. 
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PART  I  -  TURBO  PASCAL/GEM  INTERFACE 
1.  INTRODUCTION 

With  the  proliferation  of  personal  computers,  scientists  and  technologists  now  have  unprecedeiiti’d 
access  to  desktop  computing  power  at  a  level  that  just  a  few  years  ago  would  have  Iteen  avail 
able  only  from  central  mainframe  computers.  Most  notable  is  the  ready  availability  of  grapliiiv. 
procedures  within  the  high  level  languages  Turbo  Pascal  and  Turbo  C  that  considerably  eiibanii' 
the  visualisation  of  scientific  output  by  enabling  convenient  presentation  in  graphical  form.  There 
has  been,  however,  some  lag  in  the  development  of  software  suitable  for  quality  hardcopy  produc¬ 
tion.  .Although  graph  generation  software  is  available  within  various  commercial  packages  such  as 
the  GEM  software  library,  it  tends  to  be  directed  more  towards  the  business  user  whose  interest 
lies  more  in  bar  and  pie  charts  rather  than  towards  scientists  and  engineers  who  tend  to  have  a 
considerably  greater  interest  in  line  graphs. 

In  this  document,  we  describe  a  software  package  that  has  been  developed  within  Guided  Weapon- 
Division  for  interfacing  Turbo  Pascal  to  the  GEM  software  library.  .Availability  of  the  patkau'- 
allows  the  user  to  develop  and  present  graphical  data  in  standard  softcopy  form  on  an  Ili.M  l’< 
compatible  using  Turbo  Pascal,  and  at  the  same  time  to  produce  an  output  file  written  in  GKM 
format.  This  file  acts  as  the  input  to  either  GEM  Draw  or  GEM  Artline,  thus  giving  the  usm-  the 
capability  for  editing  Turbo  Pascal  generated  graphical  data  in  such  a  way  as  to  produce  liardcop> 
output  in  publication  quality  form,  incorporating  lettering  as  required  in  a  wide  variety  of  fonts  as 
well  as  whatever  artistic  touches  that  the  author  may  desire.  Even  colour  hardcopy  reprotluction 
is  possible  through  use  of  an  appropriate  colour  printer.  Experience  has  proven  that  the  TP/liK.\l 
interface  described  herein  provides  an  effective  link  between  scientific  programming/graphical  data 
analysis  and  production  of  publication  quality  graphical  output. 

The  TP/GEM  interface  and  the  essential  requirements  enabling  its  utilization  are  described  in 
Part  I  of  this  document  together  with  a  general  description  of  its  structure.  Part  II  is  wrilten 
iti  the  form  of  a  user  manual  and  includes  a  detailed  description  of  all  units  and  procedures.  .\ 
number  of  test  programs  designed  specifically  for  demonstrating  the  capabilities  of  the  interfaie 
are  included  in  the  Appendices.  A  TP/GEM  floppy  disc  is  available  on  request  to  all  interested 
users  within  DSTO.  For  best  appreciation  of  the  capabilities  of  the  interface,  it  is  recommended 
that  potential  users  request  the  disc. 


2.  HARDWARE/SOFTWARE  REQUIREMENTS 


2.1  General 

The  interface  described  in  this  document  has  been  developed  primarily  to  enable  Turbo  Ptisi  td 
programmers  to  replicate  their  softcopy  graphietd  output  in  high  quality  hardcopy  form.  The 
available  procedures,  however,  cover  rather  more  than  the  primary  graphical  output  role. 
number  of  mathematical  procedures  are  included  in  order  to  extend  the  capabilities  of  Turbo 
Pascal,  together  with  a  useful  library  of  curve  fitting  procedures.  The  interface  also  incorporates 
a  form  filling  technique  for  versatile  data  entry  at  any  point  within  a  program. 

It  is  expected  that  users  will  be  proficient  with  Turbo  Pascal  and  will  be  familiar  with  the 
graphical  design  capabilities  of  either  GEM  Draw  or  GEM  Artline.  .An  IBM  PC  compatible 
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is  required.  It  should  have  a  full  complement  of  memory  (640  KBytes)  and  have  an  EGA  or 
VG.4  Graphics  capability.  The  raw  output  is  in  the  form  of  a  GEM  format  file.  It  is  therefore 
necessary  that  sufficient  disk  storage  space  be  available  to  store  the  output  files.  The  user  must 
have  Turbo  Pascal  5  or  5.5  installed  in  order  to  use  the  graphics  procedures.  In  order  to  later 
manipulate  the  graphics  produced  in  GEM  format  the  user  must  have  GEM  Artline  or  GEM 
Draw  installed  and  a  mouse  driver.  It  is  expected  that  the  user  has  PC  or  MS-DOS  V2.1  or  a 
later  version  installed. 

2.2  Memory  Configuration 

Some  of  the  units  described  in  this  document  use  large  amounts  of  main  memory.  On  occasions. 
Turbo  Pascal  will  produce  the  message  ‘Out  of  Memory’  and  not  complete  compilation. 
number  of  techniques  can  be  used  to  resolve  this  problem  : 

1.  Check  that  you  are  compiling  to  disk  rather  than  memory  by  pressing  Alt-C  from  the  Turbo 
Pascal  Integrated  Development  Environment  (IDE).  The  pop-up  window  has  a  ‘Destination' 
option.  If  this  shows  ‘Memory’,  then  change  it  to  disk. 

2.  Reduce  the  memory  requirement.  This  may  involve  removing  all  units  that  are  not  required 
from  the  Uses  clause  in  your  program,  or  by  reducing  array  sizes. 

3.  The  Turbo  Pascal  editor  can  be  loaded  into  extended  memory  (if  installed),  saving  about 
64  KBytes  of  memory.  Consult  your  Turbo  Pascal  manual  for  details  on  how  this  may  be 
achieved. 

1.  Try  removing  memory-resident  programs  from  memory  {e.g.,  Sidekick  or  other  Terminate 
and  Stay  Resident  (TSR)  programs). 

5.  Exit  from  the  IDE  to  DOS.  Try  compiling  using  the  command  line  compiler.  For  example, 
try: 

>TPC  /L  GRAFTEST.PAS  [ENTER] 

where  /L  indicates  ‘Link  Buffer  on  Disk’,  not  memory.  For  other  TPC  options,  type; 
>TPC  [ENTER] 

All  the  command  line  options  will  appear  on  the  screen.  Alternatively,  consult  the  Turbo 
Pascal  manual. 


3.  HARD  DISK  INSTALLATION 

1.  If  you  haven’t  already  done  so,  make  a  backup  of  the  original  disk  supplied  with  this  package. 
The  easiest  way  to  do  this  is  to  use  the  DOS  diskcopy  command  (i.e.,  DISKCOPY  A:  A: 
[ENTER])  and  follow  the  prompts. 

2.  Exit  to  DOS  and  load  the  supplied  disk  into  the  floppy  disk  drive  A: 

(a)  Change  directory  to  the  level  above  where  you  want  the  files  to  be  copied  into.  For 
example,  in  order  to  place  the  files  in  a  directory  immediately  below  the  ROOT  directory. 
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type  CD\  [ENTER].  Alternatively,  if  you  prefer  to  place  them  in  a  subdirectory  of  your 
Turbo  Pascal  directory  then  type  CD\TP  (ENTER). 

(b)  Make  a  directory  into  which  the  supplied  files  can  be  copied. 

i.e..  type:  MD  PLOT  (ENTER). 

(c)  Now  move  into  your  new  directory. 

i.e.,  type:  CD  PLOT  (ENTER). 

3.  Copy  the  files  into  the  new  directory. 

i.e.,  type:  COPY  A:».»  (ENTER). 

4.  Copy  the  Turbo  Pascal  Graphics  and  Font  files  into  the  new  directory.  These  could  be  anywhere 
on  your  disk  depending  on  how  your  Turbo  Pascal  was  installed.  For  Turbo  .u,.5  iliey  an’ 
normally  located  in  the  \TP  directory.  In  this  case  type: 

COPY  \TP\  .  .BGI  [ENTER] 

COPY  \TP\  .  .CHR  [ENTER] 

Note  that  you  may  not  want  to  make  a  copy  of  these  files.  If  so,  you  can  alternatively  load  file 
PLOTLIB.P.AS  into  your  Turbo  Pascal  editor  and  locate  Procedure  ScreenSurt  (about  line 
158).  In  this  procedure  there  is  a  call  to  InitGraph  in  the  Graphics  library.  Enter  the  name  of 
the  directory  in  which  the  .BGI  and  .CHR  files  are  located  between  the  null  quotes  and  save 
the  file. 
e.g., 

InitGraphfGraphDriver.  GraphMode,  ’  ’); 
becomes 

InitGraph(GraphDriver,  Graph.Mode,  ’C:\TP\GRAPH1CS‘); 

Refer  to  the  Turbo  Pascal  reference  manual  for  details. 

5.  It  is  expected  that  users  have  already  configured  Turbo  Pascal  within  their  existing  systems. 
number  of  test  programs  (reproduced  in  the  Appendices)  are  supplied  in  order  to  illustrate  the 
capabilities  of  the  TP/GEM  interface.  It  is  suggested  that  one  of  these  files  be  loaded  into  the 
Turbo  Pascal  editor  and  compiled.  Generally  it  is  necessary  to  perform  a  BUILD  when  using 
a  unit  for  the  first  time;  i.e.,  .ALT-C,  then  B  from  the  Integrated  Development  Environment. 
Alternatively,  each  unit  can  be  picked  up  with  the  editor  and  compiled  separately. 

6.  Read  Section  4.1  before  continuing.  It  is  important  that  the  usage  of  the  type  Float, 
which  is  an  integral  part  of  the  TP/GEM  software  structure,  is  well  understood 
before  the  package  can  be  successfully  applied. 


4.  INTERFACE  STRUCTURE 

The  interface  software  is  supplied  on  a  master  diskette  in  the  form  of  a  suite  of  units,  each  of  which 
is  comprised  of  a  number  of  procedures  that  together  define  the  functional  use  of  the  unit.  The 
units  available  to  be  called  at  the  start  of  the  Turbo  Pascal  program,  all  of  which  have  a  .P.\S 
extension,  include: 
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Unit 

Role 

AXES 

Data  types  and  procedures  for  drawing  axes 

BEZIER 

Bezier  curve  generation 

COORDS 

Transformations  between  coordinate  systems 

DEFAULTS 

Sample  Default  unit 

DSM 

Discrete  smoothing  cubic  spline 

FLOATDEF 

Definition  of  the  type  Float 

FORM 

Construction  of  a  flexible  form  filling  system 

GEMMETA 

Not  intended  for  general  use 

INTERP 

Curve  data  interpolation  and  derivatives 

MATHLIB 

Supplementary  mathematiciil  library 

PLOTLIB 

Line  graph  drawing  procedures 

SBEZIER 

Smoothed  Bezier  curve  generation 

SMOOTH 

Noisy  data  smoothing 

TYPEDEF 

Global  data  types 

The  above  list  provides  a  guide  as  to  the  versatility  of  the  interface  software.  Only  the  major 
units  and  their  general  structures  are  described  here.  Part  II  provides  a  detailed  description  of 
the  variables  and  procedures  available  to  the  user.  A  number  of  examples  of  application  of  the 
procedures  for  producing  graphical  output  appear  in  the  Appendices. 

4.1  Type  Floa.t 

The  TP/GEM  interface  employs  the  type  Float  which  must  be  used  for  all  floating  point  data. 
This  is  particularly  important  when  variables  are  passed  into  a  TP/GEM  function  or  procedure, 
all  of  which  use  the  type  Float.  Float  is  defined  in  the  unit  FLOATDEF.  the  default  value  being 
Single.  This  value  may  be  modified  by  the  user  to  either  Double  or  Extended.  Examples  of 
application  of  Float  appear  throughout  the  .Appendices. 

4.2  Arrays 

Pascal  does  not  support  variable  length  arrays.  However,  Turbo  Pascal  supports  untyped  pa¬ 
rameters  that  can  be  used  to  emulate  variable  length  arrays.  When  passing  arrays  to  procedures 
that  apply  untyped  parameters  in  this  way,  the  following  rules  must  be  obeyed: 

(a)  Always  pass  the  address  of  the  first  element  of  the  data  you  wish  to  pass. 

(b)  The  array  must  have  elements  of  the  correct  type  (normally  Float). 

(c)  The  length  parameter  must  be  correct  and  less  than  the  specified  limit. 

For  example,  in  order  to  pass  the  data  contained  in  the  subarray  X[20]..X[40],  the  user  must 
specify  X[20)  as  the  array  parameter  and  21  as  the  length  parameter. 

The  limit  on  the  number  of  data  elements  passed  in  arrays  is  set  by  the  constant  MaxPt  (default 
=  1200)  in  the  unit  TYPEDEF.  This  can  be  increased  if  necessary  without  increasing  the  memory 
requirements  of  the  package. 

The  following  units  set  their  own  limits  on  the  length  of  arrays  passed  to  them; 
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INTERP  This  unit  defines  the  constant  TNArraySize  (default  =  201).  Since  the  unit  also 
declares  variables  of  the  type  ARRAY[O..TJVArraySize]  OF  Float,  altering  the 
value  of  TNArraySize  will  change  the  memory  requirements. 

SBEZIER  This  unit  defines  the  constant  InPts  (default  =  500).  Since  the  unit  also  declares 
variables  of  the  type  ARRAY[l../nPts)  OF  Float,  altering  the  value  of  InPts  will 
change  the  memory  requirements  of  the  unit. 

SMOOTH  The  unit  SMOOTH  defines  a  constant  MaxPt  (default  =  8192)  which  must  he  an 
integral  power  of  two,  and  greater  than  or  equal  to  (n  +  2  *  pts)  where  n  and  pis 
are  parameters  of  the  procedure  SmoothFt.  If  TYPEDEF. MaxPt  is  increased,  the 
user  is  advised  to  check  that  SMOOTH. MaxPt  is  still  large  enough. 


4.3  Mathematical  library 

The  common  mathematical  functions  such  as  Power.  Tan,  ArcSin,  -ArcCos  and  LoglO  that  are 
not  provided  within  Turbo  Pascal  are  included  within  the  unit  MATHLIB.  Reciprocal  conversion 
between  degree  and  radian  measure  is  also  included,  as  are  functions  for  determining  the  mantissa 
and  exponent  of  a  floating  point  number  and  a  procedure  for  locating  array  minima  and  ma.xiiiia. 


4.4  Form  filling  library 

Unit  FORM  includes  procedures  that  assist  in  the  versatile  entry  of  constant  data  into  a  program. 
It  is  a  simple  matter  to  set  up  a  form  that  is  displayed  on  the  monitor  and  which  contains  details 
of  all  the  parameters  (String,  Integer  or  Floating  Point)  that  the  user  may  wish  to  set  before 
the  program  runs.  At  this  stage  the  user  is  able  to  change  at  will  the  default  values  that  tire 
automatically  displayed,  without  the  need  to  modify  the  existing  code.  While  not  tnandutort. 
usage  of  unit  FORM  represents  a  convenience  that  becomes  quickly  apparent  when  a  program 
needs  to  be  run  many  times  with  different  starting  parameters. 

4.5  Plot  library 

The  units  that  comprise  the  plot  library  —  namely,  .AXES,  COORDS  and  PLOTLIB  —  represent 
the  heart  of  the  software  interface.  It  is  within  these  units  that  the  procedures  exist  for  creating 
the  graphical  output  in  the  form  required  by  the  user,  available  both  on  screen  and  as  a  CIF.M 
format  output  file. 

•I.5.I  Coordinate  systems 

The  Plot  Library  employs  three  different  coordinate  systems: 

World  Coordinate  Space 

This  is  the  user-defined  coordinate  system.  The  user  specifies  a  rectangular  window  withitt 
this  space  by  applying  the  DefineWorld  procedure.  Each  user-defined  window  is  calleii  ti 
World  and  is  defined  by  the  lower  left  (II)  and  upper  right  (ur)  corners.  Only  data  inside 
the  window  will  be  plotted.  The  mapping  of  a  World  into  the  plot  box  is  illustrate<l 
in  Figure  1.  The  user  can  define  many  Worlds,  thus  allowing  several  sets  of  data  to  be 
incorporated  within  the  one  plot,  as  illustrated  in  Appendix  II. 
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Figure  1  TP/GEM  coordinate  systems 


Normalized  Device  Coordinate  (NDC)  Space 

This  is  a  standard  coordinate  system  that  is  used  for  all  internal  calculations.  Each  World 
is  mapped  into  the  unit  square  as  shown  in  Figure  1.  All  curves  are  clipped  into  the  unit 
square.  All  numbering,  labelling  and  ticking  is  performed  in  NDC  space. 

Device  Coordinate  Space 

Each  output  device  is  characterized  by  its  own  coordinate  space.  The  TP/GEM  package 
supports  two  devices:  the  screen  and  the  GEM  file.  The  screen  coordinate  space  varies  with 
the  screen  type  (e.g.,  EGA  or  VGA).  The  output  plot  resides  in  a  rectangular  area.  The 
plot  box  is  centred  in  this  area,  the  relative  size  of  the  box  being  set  by  the  Scalefactor 
parameter  within  the  ScreenStart  and  GEMStart  procedures.  For  the  screen  case,  the 
rectangular  area  corresponds  to  the  entire  screen.  For  the  GEM  metafile  case,  the  area 
will  fit  onto  an  A4  page,  the  dimensions  being  set  in  the  ratio  4:3  for  the  Landscape  and 
Portrait  options  and  1:1  for  the  SquarePage  option. 

Of  prime  importance  to  users  is  the  definition  and  use  of  Worlds.  A  World  is  defined  using 
the  DeHneWorld  procedure.  The  user  specifies  the  lower  left  and  upper  right  corners  of  this 
World.  DefineWorld  returns  an  integer  value  which  represents  the  World.  This  value  should 
be  passed  to  all  the  axes  drawing,  labelling  and  plotting  procedures,  as  illustrated  within  the 
appendices.  Most  of  the  procedures  in  this  package  require  data  in  the  user-defined  World 
coordinates. 

A  few  procedures  require  data  in  NDC  coordinates  (e.g.,  PiotDashedLine).  The  user  has 
therefore  been  given  access  to  procedures  that  allow  conversion  between  the  different  coordi¬ 
nate  systems.  An  example  of  the  use  of  NDC  coordinates  is  included  in  Appendix  II.  Some 
constants  (e.g.,  TickLen)  are  defined  in  pseudo-NDC  space.  The  default  tick  length  of  0.02 
means  that  the  tick  marks  will  be  sized  at  0.02  times  the  longer  side  of  the  plot  box. 

4.6  Curve  fitting  library 

As  a  supplement  to  the  plot  library,  a  curve  fitting  library  comprised  of  the  units  BEZIER. 
SBEZIER,  DSM,  INTER?  and  SMOOTH  is  included  for  the  convenience  of  the  user,  .\s  the 
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names  suggest,  the  primary  curve  data  may  be  manipulated  before  the  output  curve  is  ploti' 
The  data  may,  for  example,  be  interpolated  or  smoothed  before  plotting.  Only  a  few  data 
points  may  be  available,  in  which  case  the  procedures  within  BEZIER  or  SBEZIER  may  In- 
applied  in  order  to  generate  smooth  connecting  curves.  The  availability  of  these  proredun  ^ 
adds  significantly  to  the  versatility  and  general  capability  of  the  TP/GEM  interface. 

4.6.1  Bezier  Curves 

GEM  Artline  supports  compound  Bezier  Curves'.  Each  segment  of  a  compound  Bezier  Ciir'.e 
is  defined  by  four  points;  two  endpoints  (iq,  J/o)  and  two  control  points  ( c, .  i;,  j 

and  {X2,y2].  The  parametric  equations  of  a  Bezier  segment  are; 


r(t)  =  I0(1  -  +  -  0^  +C2«^(1  -  0  +  X3t^ 


and 


y(t)  =  yo(l  -  <)^  +  yit(l  -  tf  +  yitUl  -  G  +  yst^ 


Two  units,  BEZIER  and  SBEZIER,  have  been  provided  for  converting  user  data  into  cominiv- 

ite  Bezier  curves.  The  advantages  of  Bezier  curves  over  normal  'curves'  composed  of  siiuiyiit 

line  segments  are; 

i)  Beziers  are  real  curves  that  are  not  composed  of  straight  line  segments. 

ii)  Compound  Bezier  Curves  are  single  entities  in  GEM  .Artline.  .'Vormal  segmented  i  iirve., 
are  limited  to  100  line  segments;  curves  of  more  than  100  segments  appear  to  (IKM 
.Artline  as  several  objects  grouped  together. 

iii)  Bezier  Curves  can  be  plotted  in  dashed  line  form.  However,  the  single  entity  nature  in 
GEM  Artline  is  then  lost. 

iv)  Bezier  curves  can  be  edited  with  ease  in  GEM  .Artline. 

Two  separate  units  have  been  provided  for  plotting  Bezier  curves; 

(1)  The  unit  BEZIER  generates  quick  curves  and  requires  less  memory  than  the  allernniivc 
unit,  SBEZIER.  BEZIER  selects  up  to  40  of  the  user-input  data  points  as  the  endpoint.^ 
of  the  Bezier  segments.  Data  points  that  are  not  chosen  as  endpoints  may  not  lie  on  the 
resulting  curve.  Furthermore,  in  some  cases  the  resulting  curve  may  not  fit  the  iniml 
data  very  well,  particularly  when  large  curvatures  are  involved. 

(2)  The  unit  SBEZIER  is  more  sophisticated.  This  unit  generates  a  smoothing  spline  for  the 
input  data  and  selects  the  endpoints  for  the  Bezier  segments  from  the  spline.  The  Bezier 
segment  equations  ate  found  by  forcing  the  first  and  second  derivatives  of  the  Bezier 
curves  to  match  those  of  the  smoothing  spUne  at  the  segment  endpoints.  The  unit 
produces  better  results  than  BEZIER  but  at  the  expense  of  greater  memory  preferem  e> 
and  slower  running  time.  Dashed  lines  of  higher  quality  (dashes  of  equal  length)  are 
generated  and  several  dashed  line  styles  are  supported. 
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PART  II  -  SOFTWARE  STRUCTURE 

In  this  second  part,  the  structure  of  the  TP/GEM  interface  is  detailed.  Sufficient  information 
is  given  to  allow  a  progammer  to  interface  between  Turbo  Pascal  and  GEM  and  to  thereby 
produce  laser  printer  output  of  the  graphical  data  displayed  initially  on  the  PC  screen. 

Examples  of  programs  that  usefully  illustrate  many  features  of  the  interface  software  are 
reproduced  within  the  Appendices  and  are  also  included  on  the  diskette.  Note  that  some 
minor  displacements  in  the  positioning  of  the  labels  appearing  in  the  raw  GEM  file  output 
have  been  corrected  within  each  of  the  graphs  shown  in  the  Appendices.  The  GEM  colours 
(which  do  not  exactly  match  the  Turbo  Pascal  softcopy  colours)  have  also  been  set  univer...iily 
to  black  (except  in  the  case  of  the  output  from  SMOOTHTEST.P.AS). 


5.  UNITS 

The  interface  structure  is  divided  into  a  number  of  units,  each  of  which  has  a  .P.\S  extension.  It  is 
advisable  to  call  the  units  FLOATDEF  and  TYPEDEF  in  ail  programs.  The  available  units  are; 


Unit 

Section 

Role 

AXES 

6.3 

Data  types  and  procedures  for  drawing  axes 

BEZIER 

6. .5 

Bezier  curve  generation 

COORDS 

0.8 

Transformations  between  coordinate  systems 

DEFAULTS 

6.12 

Sample  default  unit 

DSM 

6.7 

Discrete  smoothing  cubic  spline 

FLOATDEF 

6.10 

Definition  of  the  type  Float 

FORM 

6.1 

Construction  of  a  flexible  form  filling  system 

GEMMETA 

Not  intended  for  general  use 

INTERP 

6.9 

Curve  data  interpolation  and  derivatives 

M.ATHLIB 

6.1 

Supplementary  mathematical  library 

PLOTLIB 

6.2 

Line  graph  drawing  routines 

SBEZIER 

6.5 

Smoothed  Bezier  curve  generation 

SMOOTH 

6.6 

Noisy  data  smoothing 

TYPEDEF 

6.11 

Global  data  types 

The  complete  list  of  available  procedures  and  functions  within  the  units  is  given  in  Table  1  (in 
alphabetical  order). 

5.1  Data  types  and  global  variables 

The  units  incorporate  a  number  of  global  variables,  the  default  values  of  which  can  be  modified 
at  any  point  within  the  user  source  code,  as  illustrated  in  Appendices  I  to  111.  .Alternatively, 
the  required  changes  to  the  default  values  can  be  entered  into  the  unit  DEFAULTS. P.AS,  thus 
enabling  the  TP/GEM  interface  to  be  tailored  to  the  user  requirements.  Note  that  the  units 
in  which  the  appropriate  global  variables  are  defined  must  then  be  called  under  the  USES 
heading  within  DEFAULTS. PAS.  The  global  variables  and  data  types  are  listed  in  Tables  2  and 
3  respectively. 
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Table  1  AVAILAllI.E  PROCEDURES  AND  FUNCTIONS 


Procedure  or  Function 

I'nit 

Procedure  or  Function 

Unit 

AdvanceLine 

Form 

LogRight  Axis 

A  xes 

Arc  Cos 

MatliLib 

LogScale 

Axes 

Arc  Sin 

MathLib 

LogXAxis 

Axes 

AxisSoale 

Axes 

LogYAxis 

Axes 

ClearCoords 

Coords 

Mantissa 

MathLib 

CloseForm 

Form 

NameRightAxis 

A  xes 

CubicSplineFree 

Interp 

NameTopAxis 

A  xes 

Define  World 

Coords 

NanieXAxis 

A  xes 

Deg 

MathLib 

NameYAxis 

Axes 

Derivative 

Interp 

NDC 

Coords 

DevCoords 

<  'oords 

NDCX 

Coords 

DevX 

(  oords 

NDCY 

( 'oords 

DevY 

( 'oords 

OpenForm 

Form 

Draw  Axes 

Axes 

PlotBezier 

Bezier 

Exponent 

MathLib 

PlotCurve 

PlotLib 

FindExtrema 

MathLib 

PlotDashedLine 

PlotLib 

FitCS 

DSM 

PlotEnd 

PlotLib 

FitDS 

DSM 

PlotMarkers 

PlotLib 

FitDDS 

DSM 

Power 

MathLib 

FitDSm 

DSM 

Rad 

MathLib 

Float  Item 

1  orm 

ScreenSlart 

PlotLib 

GemWidth 

Plot  Lib 

SecondDerivative 

Interp 

GemStart 

Plot  Lib 

SetDefaults 

Defaults 

GemTextSettings 

Plot  Lib 

SmoothBezier 

S  Bezier 

Get[nteger\'alue 

Form 

SmoothFt 

Smooth 

Get  Real \'al  lie 

Form 

Stringltem 

Form 

GetStringV'alue 

Form 

Tan 

MathLib 

Integerltem 

Form 

VarLabelRight 

Axes 

LabelRightAxis 

Axes 

VarLabelTop 

Axes 

I.abelXAxis 

Axes 

World  X 

Coords 

LabelYAxls 

Axes 

World  Y 

Coords 

Log  10 

MathLib 

Table  2 

GLOBAL  VARIABLES 

\'ariable  name 

Type 

Default 

Unit 

AxisHox  = 

AxisOption; 

Full  Box 

.A  .xes 

,  .  .isColor  - 

WORD; 

Green 

.Axes 

AxisTicks  ^ 

1  it  kOption; 

TwoSldesTicked 

.A  .xes 

DeviceList  = 

OulpiitDevicesLisl; 

— 

TypeDef 

FormControl  = 

INTEGER; 

1 

Form 

LabelMargiti  = 

Float; 

0.02 

.Axes 

LogTicks  = 

LogTickOption; 

UnitsTicked 

Axes 

Side  = 

SideType; 

OutSide 

■Axes 

TickLen  = 

Float; 

0.02 

•Axes 

WSRL-GD-02/91 
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Table  3  DATA  TYPES 


Data  Type 

Imit 

Section 

.A.xisOption 

= 

(FullBox,  HalfBox); 

Axes 

6.3 

Float 

= 

REAL; 

Float  Def 

6.10 

GemFontType 

= 

(Swiss,  Swi.ssBold,  Dutch, 
DutcliBold,  Charter) 

Plot  Lib 

6.2 

GraphType 

= 

(Linear,  Log,  LogLinear, 
LincarLog); 

TypeDef 

6.11 

LineType 

= 

(Dashed,  Short  Dashed); 

PlotLib 

6.2 

LogTickOption 

= 

(DecadesTicked,  UnitsTicked); 

.Axes 

6.3 

.MarkerType 

= 

(Dot,  Plus»  Star,  Square, 

Cross.  Diamond): 

PlotLib 

6.2 

OrientType 

(Portrait,  Landscape, 
SquaroPage): 

PlotLib 

6.2 

Output  DevicesList 

Set  of  Output  DisplayTvpe; 

TypeDef 

6.11 

Out  put  Display  Type 

(Screen,  GeiuFile): 

FypeDef 

6.11 

Plot.-Vrray 

.Array[l..MaxPtl  of  Float; 

TypeDef 

6.11 

SideType 

= 

(Inside.  Outside); 

A.xes 

6.11 

TickOption 

(OneSideTickod, 

TwoSidos'l'icked ); 

Axes 

6.11 

T.N  Vector 

.Array[0..  r.N.\rr;ivSize| 
of  Float; 

DSM 

6.7 

WidthType 

(FullWidth,  llaltAVidth); 

PlotLib 

6.2 

2.  UNITS  AND  PROCEDURES 


2.1  Unit  MathLib 

The  MathLib  unit  contains  standavil  n\iUhi'malical  procedures  that  are  not  supplied  within 
Turbo  Pascal. 

Power 

Header  FUNCTIO.N  Power!;.,  b  ;  Float)  :  Float; 

Purpose  Returns  the  value  ol  n  to  the  power  of  6. 

Example  c  :=  Power(2,  x); 

c  will  assume  the  v.due  of  256  or  2®. 

Tan 


FUNCTION  Tan(x:  Float)  :  Float; 

Returns  the  tangent  of  j-  where  x  is  in  radians, 
c  ;=  Tan(tr/4); 

c  will  assume  the  value  of  1.0. 


Header 

Purpose 

Example 
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Arcsin 

Header 

Purpose 

Example 

ArcCos 

Header 

Purpose 

Example 

Rad 

Header 

Purpose 

Example 

Deg 

Header 

Purpose 

Example 

LoglO 

Header 

Purpose 

Example 

Exponent 

Header 

Purpose 

Example 

Mantissa 

Header 

Purpose 

Example 

FindExtrema 

Header 

Purpose 

Example 


FUNCTION  ArcSin(x;  Float)  :  Float; 
Returns  the  arcsine  of  x. 
c  :=  ArcSin(l/\/2); 

c  will  assume  the  value  of  )r/4  radians. 


FUNCTION  ArcCos(x;  Float)  ;  Float; 
Returns  the  arccosine  of  x. 
c  ;=  ArcCos(\/3/2); 
c  will  assume  the  value  of  jr/6  radians. 


FUNCTION  Rad(x;  Float)  ;  Float; 
Converts  degrees  to  radians, 
c  :=  Rad(45.0); 

c  will  assume  the  value  of  ir/4  radians. 


FUNCTION  Deg(x:  Float)  :  Float; 
Converts  radians  to  degrees, 
c  :=  Deg(3.14l6); 
c  will  assume  the  value  of  180°. 


FUNCTION  Logt0(x;  Float)  :  Float; 
Returns  log,o(i). 
c  ;=  LoglO(lOOO); 
c  will  assume  the  value  of  3. 


FUNCTION  Exponent(x:  Float)  ;  Float; 

Returns  the  exponent  of  a  floating  point  number, 
c  ;=  Exponent(1234..867); 

c  will  assume  the  value  of  3  (i.e.,  1234.567  =  1.234567  x  10^). 


FUNCTION  Mantissa(x;  Float)  ;  Float; 

Returns  the  mantissa  of  a  real  number, 
c  :=  Mantissa(  1234.567); 

c  will  assume  the  value  of  1.234567  (i.e.,  1234.567  =  1.234567  x  10*)- 


PROCEDURE  FindExtrema(VAR  x;  Npts  :  INTEGER; 

VAR  MinValue,  MaxValue  ;  Float;  VAR  imin,  imax  :  INTEGER); 
Finds  the  minimum  and  maximum  values  of  an  array. 

Also  returns  the  index  positions  where  the  minima  and  ma.xima  were  foiiii 
FindExtrema(PlotData,  100,  Minimum,  Maximum,  MinPos,  MaxPos); 
Returns  the  maximum  and  minimum  values  of  the  first  100  elements 
of  the  array  PlotData.  The  minimum  and  maximum  values  are 
located  at  index  positions  MinPos  and  MaxPos  respectively. 
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6.2  Unit  PlotLib 

Unit  PlotLib  conttiins  line  graph  drawing  procedures  to  screen  and/or  to  a  GEM  metafile. 


6.2.1  Enumerated  data  types 


OrientType 

GemFontType 

LineType 

VVidthType 

MarkerType 


(Portrait,  LandScape,  SquarePage); 

(Swiss,  SwissBold,  Dutch,  DutchBold,  Charter): 
(Dashed,  ShortDashed); 

(FullWidth,  HalfWidth); 

(Dot,  Plus,  Star,  Square,  Cross,  Diamond); 


6.2.2  Mandatory  Calls 

The  following  list  includes  the  mandatory 

ScreenStart/GemStart 
Draw  Axes 
DefineWorld 

LabelXAxis  or  LogXAxis(unit  Axes) 
LabelYAxis  or  LogYAxis(unit  Axes) 
PlotCurve  (or  PlotBezier) 

Plot  End 


calls  for  drawing  a  curve: 

:  Initialise  screen  or  GEM  file 

Draws  rectangle  around  graph  area 
:  Define  World  (user)  coordinates 

Draw  tick  marks  and  number  X  .Axis 
Draw  tick  marks  and  number  Y  .Axis 
Draw  a  curve 
Finish  plot 


6.2..3  Available  procedures  and  functions 

GemWidth  Sets  line  width  for  Gem  polylines 

GemStart  Initialises  Gem  plot  file  (GEMfile) 

GemTextSettings  Sets  font  size  and  type  for  Gem  file  text 

PlotCurve  Plots  the  contents  of  an  array 

Plot  Dashed  Line  Sets  the  color,  dash  type  and  line  width 

PlotEnd  Called  when  plot  is  complete 

PlotMarkers  :  Used  to  draw  data  points 

ScreenStart  :  Used  to  initialise  and  scale  plotting  to  screen 


6.2.4  Procedure  and  function  headers 

Reference  should  be  made  to  the  above  enumerated  data  types. 

PROCEDURE  GemWidth(LineWidth  :  INTEGER); 

LineWidth  Width  of  GEM  line  (default  value  =  50  (0.5  mm). 


PROCEDURE  GemStart(FileName  :  STRING;  Orientation  :  OrientType: 
ScaJeFactor  :  Float); 

FileName  Name  of  GEM  output  file.  Call  this  file  from  GEM  draw 

or  GEM  Artline  to  access  your  graphics  output. 
Orientation  Orientation  of  output  as  per  type  OrientType. 

ScaleFactor  Defines  the  fraction  of  the  available  plot  dimension 

within  the  plot  rectangle. 
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PROCEDURE  GemTextSettings  (Font  ;  GcmFontType;  FontSize  :  INTEGER); 

Font  :  Specifies  font  type  as  per  typeGemFontType  (default  =  Swiss). 

FontSize  :  Point  size  for  font  (default  =  17  pt). 

PROCEDURE  PlotCurve(Color  ;  WORD;  NPoints  ;  INTEGER;  GtaphKind  ;  GtaplvType. 
VAR  X,  y;  World  :  INTEGER); 

Color  Color  of  plotted  curve  (standard  Turbo  Pascal  screen  colors). 

NPoints  Number  of  points  in  curve. 

GraphKind  (Linear,  Linearlog,  LogLinear,  Log). 

X,  y  Arrays  containing  the  curve  (i,  y)  coordinates. 

World  :  The  World  associated  with  the  curve. 

PROCEDURE  PlotDaahedLine( Color  :  WORD;  LineForm  :  LineType; 

DashedVVidth  :  WidthType;  xl,  yl,  x2,  y2  :  Float); 

Color  Color  of  plotted  curve. 

LineForm  Dashed  or  ShortDashed. 

DashedWidth  FullWidth  or  HalfWidth. 

xl,  yl  One  end  point  of  the  line  (NDC  Coordinates). 

x2,  y2  Other  end  point  of  the  line  (NDC  Coordinates). 

PROCEDURE  PlotEnd; 

PROCEDURE  PIotMarkers(Color  :  WORD;  Kind  :  MarlcerType; 

Size  ;  Float;  NPts  •.  INTEGER;  GtaphKind  GraphType;  Vat  x.  y;  World  ;  lNTEGF.il ): 

Color  Color  of  plotted  curve. 

Kind  Type  of  marker  (see  Marker  type  above). 

Size  :  Size  of  marker  relative  to  standard  (  =  1)  device  size. 

Npts  ;  Number  of  markers  (or  points)  to  plot. 

GraphKind  :  (Linear,  Linearlog,  LogLinear,  Log). 

x,  y  Arrays  containing  the  x  and  y  values  of  the  points  to  plot. 

World  The  World  associated  with  the  curve. 

PROCEDURE  ScreenStart(ScaleFactor  :  Float); 

ScaleFactor  Fraction  of  full  screen  dimensions  allocated  for  the  plot  rectangle. 


6.3  Unit  Axes 

Unit  Axes  contains  data  types  and  procedures  for  drawing  axes  on  graphs.  The  unit  also  contains 
an  autoscaling  procedure  that  may  be  applied  separately  to  all  curve  axes. 

6.3.1  Available  structures 

6.3. 1. 1  Data  types 

The  following  descriptive  datatypes  are  used  to  describe  the  axes  required: 
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TickOption 
SideType 
LogTickOption 
AxisOption 
Array  100 


(OneSideTicked,  TwoSidesTicked); 
(Inside,  Outside); 

(DecadesTicked,  UnitsTicked); 
(FullBox,  HalfBox); 
ARRAY[1..100)  OF  Float; 


6.3. 1.2  Global  variables  and  defaults 
AxisOption; 
WORD; 


AxisBox 

AxisColor 

AxisTicks 

LabelMargin 

LogTicks 

Side 

TickLen 


TickOption; 

Float; 

LogTickOption; 

SideType; 

Float; 


FullBox 

Green 

OneSideTicked 

0.02 

UnitsTicked 

Outside 

0.02 


Note  that  a  grid  can  be  established  by  setting  Side  to  InSide  and  TickLen  to  1,0.  The 
width  of  the  grid  lines  may  be  set  through  use  of  the  procedure  GEMWidtb,  as  illustrated 
within  Appendix  III. 


6.3.2  Available  procedures 


AxisScale 

LabelRightAxis 

LabelXAxis 

LabelYAxis 

LogRightAxis 

LogSeale 

LogXAxis 

LogYAxis 

NameRightAxis 

NameTopAxis 

NameXAxis 

•NameYAxis 

VarLabelRight 

VarLabelTop 


Provides  autoscaled  values  from  specified  min  and  max  values. 
Draws  ticks  and  numbers  on  the  right  axis. 

Draws  ticks  and  numbers  on  the  bottom  axis  and 
optionally  draws  ticks  on  the  top  axis. 

Draws  ticks  and  numbers  on  the  left  axis  and 
optionally  draws  ticks  on  the  right  axis. 

As  for  LabelRightAxis  except  with  log  scale. 

Provides  autoscaled  logarithmic  values  from  specified 
maximum  value  and  required  number  of  decades. 

As  for  LabelXAxis  except  with  log  scale. 

As  for  LabelYAxis  except  with  log  scale. 

Writes  a  name  on  the  right  axis. 

Writes  a  name  on  the  top  axis. 

Writes  a  name  on  the  bottom  axis. 

Writes  a  name  on  the  left  axis. 

Draws  user-defined  ticks  and  labels  on  the  right  axis. 

Draws  user-defined  ticks  and  labels  on  the  top  axis. 


6.3.3  Procedure  headers 


PROCEDURE  AxisScale(AxisnameStr  ;  STRING;  VAR  steps  :  INTEGER; 
VAR  min,  max  :  Float;  VAR  Decimals  ;  INTEGER); 


AxisnameStr 
Steps 
min,  max 

Decimals 


Identifying  name  for  axis. 

Number  of  steps  required  on  the  axis. 

Minimum  and  maximum  values  (input  as  specified, 
output  as  autoscaled  or  user-selected). 

Number  of  displayed  decimals. 
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PROCEDURE  LabelRightAxis(LabelColor  ;  Word;  Ylntcrvals,  decimals  :  INTEGER; 
World  :  INTEGER); 

LabelColor  Color  assigned  to  the  ticks  and  numbers. 

YIntervals  Total  intervals  between  right  hand  axis  ticks. 

Decimals  Number  of  displayed  decimals. 

World  The  world  associated  with  the  plotted  curve. 

PROCEDURE  LabelXAxis( LabelColor  ;  WORD;  Xlntervals,  Decimals  :  INTEGER: 
World  :  INTEGER); 

LabelColor  Color  assigned  to  the  ticks  and  numbers. 

Xlntervals  Total  intervals  between  bottom  axis  ticks. 

Decimals  Number  of  displayed  decimals. 

World  The  world  associated  with  the  plotted  curve. 

PROCEDURE  LabelYAxis( LabelColor  :  WORD;  YIntervals,  Decimals  :  I.N’TEGER: 
World  ;  INTEGER); 

LabelColor  Color  assigned  to  the  ticks  and  numbers. 

YIntervals  Total  intervals  between  left  hand  axis  ticks. 

Decimals  Number  of  displayed  decimals. 

World  The  world  associated  with  the  plotted  curve. 

PROCEDURE  LogRightAxis( LabelColor  :  WORD;  World  ;  INTEGER): 

LabelColor  Color  assigned  to  the  right  hand  ticks  and  numbers. 

World  The  world  associated  with  the  plotted  curve. 

PROCEDURE  LogScale(  AxisnameStr  ;  STRING;  VAR  min.  max  :  Float: 

NoDecades:  INTEGER); 

AxisnameStr  Identifying  name  for  axis, 

min,  max  Minimum  and  maximum  values 

—  max  is  input  as  specified, 

output  as  autoscaled  or  user-selected. 

—  min  is  calculated  from  NoDecades  and  ma.v 

NoDecades  :  Number  of  decades  required  on  log  axis. 

PROCEDURE  LogXAxisf LabelColor  :  WORD;  World  :  INTEGER); 

LabelColor  Color  assigned  to  the  bottom  axis  ticks  and  numbers. 

World  The  world  associated  with  the  plotted  curve. 

PROCEDURE  Log YAxis( LabelColor  :  WORD;  World  :  INTEGER); 

LabelColor  Color  assigned  to  the  left  hand  axis  ticks  and  numbers. 

World  The  world  associated  with  the  plotted  curve. 

PROCEDURE  NameRightAxi8( Color  :  WORD;  YNameStr  ;  STRING), 

Color  Color  assigned  to  the  right  axis  name. 

YNameStr  Right  axis  name. 
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PROCEDURE  NameTopAxi«( Color  :  WORD;  XNameStr  :  STRING): 

Color  Color  assigned  to  the  top  axis  name. 

XNameStr  Top  axis  name. 

PROCEDURE  NameXAxia(Color  :  WORD;  XNameStr  :  STRING); 

Color  Color  assigned  to  the  bottom  axis  name. 

XNameStr  Bottom  axis  name. 

PROCEDURE  NameYAx!s(Color  :  WORD;  YNameStr  :  STRING): 

Color  :  Color  assigned  to  the  left  hand  axis  name. 

YNameStr  :  Left  hand  axis  name. 

The  following  two  procedures  give  the  opportunity  for  user-defined  labelling.  There  can  be 
cases  where  an  axis  needs  to  be  labelled  linearly  in  terms  of  a  primary  parameter  on  one  .side, 
and  nonlinearly  in  terms  of  a  functionally  related  parameter  on  the  other  side,  as  illustrated  in 
.Appendix  II.  The  user  must  specify  two  arrays;  MyLabels,  for  defining  the  required  labels  to 
be  displayed  nonlinearly  on  the  axis  of  interest,  and  TickPosn,  which  is  derived  by  converting 
the  MyL&beh  values  into  the  corresponding  linear  parameter  values,  in  accordance  with  the 
required  nonlinear  functional  form.  For  example,  if  the  primary  parameter  is  wavelength  A 
to  be  displayed  linearly  on  the  lower  x-axis  as  in  the  example  of  .Appendix  11.  and  the  user 
wishes  to  display  the  energy 


1.240 

A 


(!) 


across  the  upper  x-axis,  then  the  array  MyLabels  defining  the  required  gradations  of  E 
is  formed,  followed  by  the  TickPosn  array  which  defines  the  complementary  values  of  A 
obtained  from  the  inverse  form  of  the  functional  expression  (1)  above.  This  example  is  coded 
in  Appendix  II  where  it  may  be  noted  that  both  TickPosn  and  AfyLabe/s  are  specified  as 
arrays  of  Float. 

PROCEDURE  VttrLabelRight(LabelColor  ;  WORD;  NTicks  .  INTEGER. 

VAR  TickPosn,  MyLabels;  Decimals  :  INTEGER;  YMin,  YMax  ;  Float); 

LabelColor  Color  of  right  axis  label,  ticks  and  numbers. 

NTicks  Number  of  ticks  to  display  on  right  axis. 

TickPosn  Array  of  Y  coordinates  where  ticks  are  to  be  placed. 

MyLabels  Array  of  labels  for  tick  marks. 

Decimals  Number  of  displayed  decimals. 

YMin,  YMax  Minimum  and  Maximum  values  of  Y  for  the  right  axis. 

PROCEDURE  VarLabelTopi LabelColor  .  WORD;  NTicks  ;  INTEGER: 

VAR  TickPosn,  MyLabels;  Decimals  ;  INTEGER;  XMin,  XMax  :  Float); 

LabelColor  Color  of  right  axis  label,  ticks  and  numbers. 

NTicks  Number  of  ticks  to  display  on  right  axis. 

TickPosn  Array  of  Y  coordinates  where  ticks  are  to  be  placed. 

MyLabels  Array  of  labels  for  tick  marks. 

Decimals  Number  of  displayed  decimals. 

XMin,  XMax  Minimum  and  Maximum  values  of  X  for  the  top  axis. 
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6.4  Unit  Form 

Unit  Form  provides  procedures  for  flexible  data  entry  on  the  same  basis  as  filling  out  a  form.  Tlio 
user  constructs  the  desired  layout  of  the  form  by  application  of  the  available  procedures.  The 
input  data  may  be  in  the  form  of  strings,  integers  or  floating  point  numbers.  Use  of  the  CURSOR 
or  ENTER  keys  allows  movement  between  different  fields  and  pves  the  opportunity  for  the  user 
to  change  the  programmed  default  data.  Data  entered  in  incorrect  numeric  format  is  ignored. 
The  data  on  the  form  is  accepted  when  either  of  the  PGDN  or  ESC  keys  is  pressed.  Refer  to 
Appendix  VII  for  a  detailed  description  of  unit  Form  within  the  test  program  FORMTEST.P.IS. 

6.4.1  Available  procedures 

AdvanceLine  Inserts  a  blank  line  in  the  form. 

CloseForm  Indicates  that  the  layout  is  complete. 

Floatltem  ;  Places  a  field  for  floating  point  input  at  the  required  location. 

GetIntegerValue  Inputs  an  integer. 

GetRealValue  Inputs  a  floating  point  number. 

GetStringValue  Inputs  a  string. 

Integerltem  Places  a  field  for  an  integer  input  at  the  required  location. 

OpenForm  Starts  construction  of  the  layout. 

Stringitem  Places  a  field  for  a  string  input  at  the  required  location. 

6.4.2  Procedure  headers 
PROCEDURE  AdvanceLine; 

PROCEDURE  CloseForm; 

PROCEDURE  Float Item( Float Val  ;  Float;  Decimals  :  INTEGER;  ItemSir  ;  STRING); 

Float Val  Float  variable  name  for  the  field. 

Decimals  Number  of  displayed  decimals. 

ItemStr  :  Prompt  associated  with  the  variable  name. 

PROCEDURE  GetIntegerValue  (VAR  IntegerValue  :  I.NTEGER); 

IntegerValue  Returns  the  integer  value  input  at  the  keyboard. 

PROCEDURE  GetRealValue  (VAR  RealVaiue  :  Float); 

RealValue  Returns  the  real  (or  Float)  value  input  at  the  keyboard. 

PROCEDURE  GetStringValue  (VAR  StringVal  ;  STRING); 

StringVal  Returns  the  string  input  at  the  keyboard. 

PROCEDURE  Integerltem  (IntegerVal,  Width  ;  INTEGER;  ItemStr  ;  STRING); 

IntegerVal  Integer  variable  to  be  used  with  this  field. 

Width  Width  of  displayed  field. 

ItemStr  Prompt  associated  with  the  variable  name. 
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PROCEDURE  OpenForiii(HeadingStr  :  STRING); 

HeadingStr  Title  for  top  of  form. 

PROCEDURE  Stringltem  (StringVal  :  STRING;  width  :  INTEGER;  ItemStr  :  STRING); 

StringVal  String  variable  name  to  be  used  with  this  field. 

ItemStr  Prompt  to  associate  with  this  variable  name. 

6.5  Units  Bezier  and  SBezier 

These  units  contain  procedures  for  generating  Bezier  curves.  Bezier  curves,  which  are  supported 
by  GEM  Artline  but  not  by  GEM  Draw,  are  used  for  drawing  smooth  curves  through  a  limited 
number  of  data  points.  An  example  of  use  of  PlotBezier  can  be  found  in  file  TESTDSM.P.^S 
(reproduced  within  Appendix  V).  An  example  of  use  of  SmoothBezier  is  to  be  found  in  file 
GRAFTEST.PAS  (reproduced  within  Appendix  III). 


6.5.1  Data  type  (Available  only  with  SBEZIER) 

DashedLineComponent  =  (LongDash.  ShortDash.  Dot,  LongSpace.  ShortSpace): 


6.5.2  Available  procedures 

PlotBezier  (Unit  Bezier)  Plots  a  Bezier  curve. 

SmoothBezier  :  (Unit  SBezier)  Fits  to  a  smoothing  Bezier  curve. 

SBezier  provides  a  better  fit,  the  lengths  of  the  dashed  segments 
tend  to  be  more  uniform,  but  it  is  slower  and  uses  more  memory. 


6.5.3  Procedure  headers 


PROCEDURE  PlotBezier  (Color,  NPts,  NDashes  :  INTEGER;  GraphKind  :  GraphType; 
VAR  X,  y,  wx,  wy;  World  :  INTEGER); 

Color  Color  of  the  Bezier  curve  plot. 

NPts  Number  of  input  data  points. 

NDashes  Number  of  dashes  if  dashed  curve  is  required 

(=  0  for  a  continuous  curve). 

GraphKind  (Linear,  Linearlog,  LogLinear,  Log). 

X,  y  Arrays  containing  the  (x,y)  input  curve  data, 

wx,  wy  Work  arrays  (must  be  of  dimension  NPts). 

World  The  world  associated  with  the  Bezier  curve. 


PROCEDURE  SinoothBezier(Color  :  WORD;  NPts,  LineType  :  INTEGER; 
GraphKind  :  GraphType;  VAR  x,  y;  World  :  INTEGER); 
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Color 

NPts 

LineType 


GraphKind 
X,  y 


Color  of  the  Bezier  curve  plot. 

Number  of  input  data  points. 

1  —  Long  dashes  and  short  spaces. 

2  —  Short  dashes  and  short  spaces. 

3  —  Dots  with  short  spaces. 

4  —  Medium  dashes  with  short  spaces. 

5  —  Long  dash,  short  space,  dot,  short  space. 

6  —  Short  dashes  with  long  spaces. 

7  —  Dots  with  long  spaces. 

(Linear,  Linearlog,  LogLinear,  Log). 

Arrays  containing  the  (x,y)  input  curve  data. 


6.6  Unit  Smooth 


This  unit  is  available  for  smoothing  large  amounts  of  noisy  curve  data.  The  data  values  in  tlm 
X  array  must  be  equally  spaced.  The  procedure  SmoothFt^  operates  by  applying  a  Fast  Foiuici 
Transform  to  the  data  in  the  y  array,  followed  by  a  low  pass  filter.  Linear  trends  are  retained-. 

Program  SMOOTHTEST.PAS  (Appendix  VI)  contains  an  example  of  application.  In  this  pro¬ 
gram  a  sine  curve  is  plotted,  random  noise  is  added,  and  the  original  curve  compared  to  the 
smoothed  curve  obtained  by  application  of  SmoothFt. 


6.6.1  Procedure  header 

PROCEDURE  SmoothFt(VAR  y;  n  :  INTEGER;  Pts  :  Float): 

y  ;  Unsmoothed  input  and  smoothed  output  y  array, 

n  :  Number  of  elements  in  array  y  (maximum  =  8192). 

Pts  Number  of  points  over  which  the  data  is  to  be  smoothed 

0  —  No  Smoothing 

>  0.5n  —  Both  data  and  noise  are  smoothed 


6.7  Unit  DSM 

Unit  DSM,  which  is  useful  for  smoothing  moderate  amounts  of  noisy  curve  data,  is  based  on 
the  algorithm  ALG547  —  Discrete  Smoothing  Cubic  Spline.  Credits  and  references  are  included 
in  the  Unit  source  code.  An  example  of  use  can  be  found  in  file  TESTDSM.P.4S  (reproduced 
within  Appendix  V). 

Within  the  DSM  procedures,  parameters  NInterp,  xs,  ys.  b,  c  and  d  define  the  spline.  The 
procedures  should  be  called  in  the  following  order: 

1.  FitDSm  (To  evaluate  the  spline) 

2.  FitCS  (optional) 

3.  FitDS  (optional) 

4.  FitDDS  (optional) 
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6.T.1  Available  procedures 

PitCS  Evaluates  a  cubic  spline  at  each  point  in  an  array. 

FitDS  Evaluates  the  first  derivative  of  the  cubic  spline. 

FitDDS  Evaluates  the  second  derivative  of  the  cubic  spline. 

FitDSm  Computes  the  discrete  natural  cubic  spline  defined  over  an  interval. 


6.7.2  Procedure  headers 


PROCEDURE  Fites  (NPts,  NInterp  :  INTEGER;  VAR  x,  y,  xs,  ys,  b,  c,  d): 


NPts 

NInterp 

X 

y 

xs 

ys 

b 

c 

d 


Number  of  points  to  interpolate. 

Number  of  nodes  (x)  and  data  values  (y). 

Array  of  x  coordinates  for  interpolation. 

Returned  values  of  y. 

.4rray  of  floats  containing  nodes  (i(i]  <  i[i  +  1]). 

Array  of  floats  containing  the  smoothed  values  of  data  j/[i]. 
Array  of  floats  containing  coefficients  for  terms  {t  -  j:[i]) 
(where  I  is  the  interval  parameter  defined  within  ALG.o47). 
Array  of  floats  containing  coefficients  for  terms  (t  -  x(!])". 
Array  of  floats  containing  coefficients  for  terms  (f  -  x[!]p. 


PROCEDURE  FitDS  (NPts,  NInterp  .  INTEGER;  VAR  x.  y.  xs.  ys,  b.  c.  d); 


NPts 

.NInterp 

X 

y 

XS 

ys 

b 

c 

d 


Number  of  points  to  interpolate. 

Number  of  nodes  (x)  and  data  values  (y). 

.^rray  of  x  coordinates  for  interpolation. 

Returned  values  of  dyjdx. 

Array  of  floats  containing  nodes  (i[i)  <  i[i  +  1]). 

Array  of  floats  containing  the  smoothed  values  of  data  j/[i]. 
Array  of  floats  containing  coefficients  for  terms  ((  -  xfi]) 
(where  t  is  the  interval  parameter  defined  within  .4LG547). 
Array  of  floats  containing  coefficients  for  terms  ((  -  x[i))'. 
Array  of  floats  containing  coefficients  for  terms  ((  -  x[i])^. 


PROCEDURE  FitDDS  (NPts,  NInterp  :  INTEGER;  VAR  x,  y,  xs,  ys,  b,  c,  d); 


NPts 

NInterp 

X 

y 

xs 

ys 

b 

c 

d 


Number  of  points  to  interpolate. 

Number  of  nodes  (x)  and  data  values  (y). 

Array  of  x  coordinates  for  interpolation. 

Returned  values  of  d'^y/dx^. 

Array  of  floats  containing  nodes  (i[i]  <  i[i  +  1]). 

Array  of  floats  containing  the  smoothed  values  of  data  j/[i]. 
Array  of  floats  containing  coefficients  for  terms  ((  -  x[i]) 
(where  t  is  the  interval  parameter  defined  within  ALG547). 
Array  of  floats  containing  coefficients  for  terms  (t  -  i[i])^. 
Array  of  floats  containing  coefficients  for  terms  (f  -  i[t])^. 
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PROCEDURE  FUDSm  (h,  Rho  :  Float;  n  ;  INTEGER;  VAR  x,  y,  Wgs,  ys.  b,  c.  d); 


h 

Rho 


n 

X.  y 

Wgs 

ys 

b 


c 

d 


Step  size  used  for  discrete  cubic  spline. 

Positive  parameter  for  varying  smoothness  of  fit. 

—  Large  Rho  emphasizes  smoothness 

—  Small  Rho  emphasizes  fitting 
Number  of  nodes  (x)  and  data  values  (y). 

Array  of  floats  values  containing  nodes  (i[j]  <  i[i  +  1]). 

Array  of  floats  containing  weights  W5s[i]  corresponding  to  the 
data  (i[i],  j([t]). 

Array  of  floats  containing  the  smoothed  values  of  data  y[>J. 
Array  of  floats  containing  coefficients  for  terms  (t  —  i[i]) 
(where  t  is  the  interval  parameter  defined  within  ALG.547). 
Array  of  floats  containing  coefficients  for  terms  {t  - 
.Array  of  floats  containing  coefficients  for  terms  ((  - 


6.8  Unit  Coords 

This  unit  contains  routines  for  performing  transformations  between  coordinate  systems,  .A 
maximum  of  10  worlds  and  10  device  spaces  can  be  used.  Examples  of  use  of  iho  unit  can  be 
found  throughout  the  example  files  given  in  the  appendices  to  this  document. 

6.8.1  Available  procedures 

ClearCoords  :  Deletes  all  world  and  device  spaces. 

DefineWorld  Defines  a  world  coordinate  system.  V' odd  points  can 

be  mapped  into  N  DC  space. 

DevCoords  Converts  a  point  in  NDC  space  into  device  space. 

DevX  ;  Converts  an  x  coordinate  in  NDC  space  into  device  space. 

DevY  Converts  a  y  coordinate  in  NDC  space  into  device  space. 

NDC  Converts  a  point  in  world  space  into  NDC  space. 

NDCX  :  Converts  an  x  coordinate  in  World  space  into  NDC  space 

NDCY  :  Converts  a  y  coordinate  in  World  space  into  NDC  space 

WorldX  :  Converts  an  x  coordinate  in  NDC  space  into  world  space. 

WorldY  :  Converts  a  y  coordinate  in  NDC  space  into  world  space. 


6.8.2  Procedure  and  function  headers 
PROCEDURE  ClearCoords; 

PROCEDURE  DefineWorld(VAR  World  ;  INTEGER;  llx.  lly,  urx,  ury  ;  Float); 

World  Index  name  for  new  world. 

llx,  lly  Lower  left  coordinate  of  world  coordinate  space. 

urx,  ury  Upper  right  coordinate  of  world  coordinate  space. 

PROCEDURE  DevCoords( Device  :  INTEGER;  u,  v  :  Float;  VAR  x,  y  .  INTEGER). 


Device 

u,v 

X.  y 


Index  number  for  device. 

NDC  device  coordinates. 

NDC  coordinate  converted  into  device  coordinates. 
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PROC  EDURE  DevX( Device  :  INTEGER;  u  :  Float)  :  INTEGER; 

Device  Index  number  for  device. 

n  ;  X  coordinate  in  NDC  space. 

—  returns  x  coordinate  converted  to  device  space. 

PROCEDURE  DevY(Device  :  INTEGER;  u  :  Float)  :  INTEGER; 

Device  Index  number  for  device. 

11  :  y  coordinate  in  NDC  space. 

—  returns  y  coordinate  converted  to  device  space. 

PROCEDURE  NDC(World  :  INTEGER;  x,  y  :  Float;  VAR  u,  v  ;  Float); 

World  ;  Index  number  for  world. 

-X.  y  :  Point  in  world  coordinate  space. 

11,  V  :  Point  (x,y)  converted  to  NDC  coordinates. 

FUNCTION  NDCXfWorld  :  INTEGER;  x  :  Float)  :  Float; 

World  Index  number  for  world. 

X  X  value  in  world  space. 

—  returns  x  converted  to  NDC  space. 

FUNCTION  NDCYfWorld  :  INTEGER;  y  ;  Float)  ;  Float; 

World  Index  number  for  world, 

y  ;  y  value  in  world  space. 

—  returns  y  converted  to  NDC  space. 

FUNCTION  WorldX( World  ;  INTEGER;  x  ;  Float)  ;  Float; 

World  Index  number  for  world. 

X  X  coordinate  in  NDC  space. 

—  returns  x  converted  to  world  space. 

6.9  Unit  Interp 

This  unit  provides  procedures^  for  returning  the  first  and  second  derivatives  of  a  cubic  spline. 
It  also  provides  a  routine  for  constructing  a  curve  through  a  given  set  of  data  points  (where  the 
second  derivative  at  the  endpoints  is  assumed  to  be  nero).  The  listing  within  the  lest  program 
TE.Sl  .M.ARKER.PAS  at  Appendix  V  provides  an  example  of  the  use  of  this  unit. 

fi.9.1  Data  types  and  constants 

CONST  ;  TNArraySize  =  201;  (limits  array  size). 

TYPE  TNVector  =  .ARR,AY[0.. TNArraySize]  of  Float; 


6.9.2  Available  procedures 


Derivative 

SecondDerivative 

CubicSplineFree 


Returns  the  first  derivative  of  a  spline. 

Returns  the  second  derivative  of  a  spline. 

Constructs  a  smooth  curve  through  a  given  set  of  data  points. 
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6.9.3  Procedure  headers 


PROCEDURE  Derivative  (N'limPoints  :  INTECJER:  VAR  XData;  Nunilnter  :  INTECEli: 
\'AR  Xliiter,  YInter); 


N  urn  Points 

XData 

Numinter 

Xlnter 

YInter 


Number  of  data  points. 

Summed  Cubic  Spline  Free. 

Number  of  points  where  the  derivative  is  required. 
X  values  where  the  derivative  is  required. 
Derivative  at  Xi  (i  =  E.Numlter). 


PROCEDURE  SecondDerivative  (NumPoints  :  INTEGER:  V.AR  XData: 
Numinter  ;  INTEGER:  VAR  Xlnter,  YInter); 


NumPoints 

XData 

Numinter 

Xlnter 

YInter 


Number  of  data  points. 

Summed  Cubic  Spline  Free. 

.Number  of  points  where  the  derivative  is  required. 
X  values  where  the  derivative  is  required. 

Second  derivative  at  x,  (i  =  1  ...Numlter). 


PROCEDURE  CubicSplineFVee  (NumPoints  ;  INTEGER:  V.AR  .XData.  YDaia: 
Numinter  :  I.NTEGER:  V.AR  Xlnter,  YInter:  V.AR  Error  :  BA'TE); 


.NumPoints 

XData 

YData 

.Numinter 

.Xlnter 

YInter 

Error 


.Number  of  data  points. 

Input  X  data  values. 

Input  y  data  values. 

•Number  of  interpolations. 

X  coordinates  of  points  at  which  to  interpolate. 
Interpolated  values  at  Xlnter. 

0  —  .No  Error. 

1  —  X  values  of  data  points  not  unique. 

'2  —  X  values  of  the  data  points  not  in  a.srpnding  order. 
3  —  NumPoints  <  2. 


6.10  Unit  FloatOef 


The  purpose  of  this  unit  is  to  allow  the  user  to  specify  the  type  Float,  either  as  SINGI.K. 
DOUBLE  or  EXTENDED.  The  default  value  is  SINGLE.  The  unit  must  be  included  with  all 
Turbo  Pascal  code  that  makes  use  of  the  TP/GEM  .software  interface  and  is  edited  as  required 
by  the  user. 


8.11  Unit  TypeDef 


This  unit  should  be  included  with  all  source  programs  that  incorporate  any  code  from  the 
TP/GEM  interface.  There  are  no  functions  or  procedures,  only  data  types  and  variables. 


6. II. I  Data  types 
PlotArray 

OutputDisplayTypc 

OutputDevicesList 

GemColorMap 

GraphType 


ARRAY[l..MaxPt]  OF  Float; 
(Screen,  GemFile); 

•Set  of  OutputDisplayTypc; 
ARRAY[0..1.5]  OF  INTEGER; 
(Linear,  Log,  LogLinear.  LinearLog); 
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CONST  GemColors  :  GemColorMap  =  (0,12,11,14,10,15,13.8,9,4,3,5,2,7.6,1); 


6.11.2  Global  variables 


DeviceList 

ErrorStr 

ScreenDevice 

Full  Screen  Device 

GEMDevice 

GEMPage 


OutputDevicesList; 

STRING; 

INTEGER; 


—  Screen  device  space  covered  by  the  specified  plot  rectangle 
INTEGER; 

—  Device  space  covered  by  full  screen 
INTEGER; 

—  GEM  device  space  covered  by  the  specified  plot  rectangle 
INTEGER; 

—  Device  space  covered  by  full  GEM  plot  area 


The  user  has  the  option  of  writing  data  to  the  PC  screen  but  not  to  the  GE.VI  output  file. 
The  latter  option  is  often  more  conveniently  accomplished  during  manipulation  within  GE.M 
Draw  or  GEM  Artline.  The  ScreenDevice  variable  is  useful  as  an  aid  for  positioning  the 
required  data  on  the  screen.  For  example: 

KoveTo(OevX(ScreenDevic«  .O.Ol'i'lO  ,DevY(ScreeiU)evica,l .  0)+20)  ; 

0utText( 'Filter  Cutoff  :  ’); 

Str (FilterOff :3 : 1 , VritaStr) ; 

QutTeit(WriteStr) ; 

QutTextC  microns’); 


will  write  the  required  information  near  to  the  top  left  hand  corner  of  the  specified  plot 
rectangle.  A  number  of  examples  of  the  application  of  ScreenDevice  appear  in  .Appendix  II. 

6.12  Unit  Defaults 

The  purpose  of  this  unit  is  to  allow  the  user  to  tailor  the  TP/GEM  interface  to  suit  his/her  own 
personal  preferences.  .Any  value  of  a  global  variable  written  into  the  procedure  SetDefau/ts  will 
overwrite  the  default  value.  Note  that  the  unit  in  which  the  variable  is  defined  must  be  called 
under  the  USES  heading  within  DEFAULTS. PAS  which,  in  turn,  must  be  called  within  the  u.sor 
source  program.  Note  further  that  the  values  of  the  global  variables  may  be  changed  at  any 
point  within  the  user  source  program,  as  illustrated  for  example  within  Appendices  I  to  III. 


PROCEDURE  SetDefaults; 


2S 
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Appendix  I 

Test  Program:  EASYPLT.PAS 


This  is  a  simple  program  designed  to  illustrate  the  basic  use  of  the  TP/GEM  interface,  with  data  im¬ 
ported  from  an  external  data  file.  Note  the  example  of  the  change  of  the  global  variable  AxisColor. 
This  user  preference  could  optionally  have  been  stored  in  the  user  defaults  unit  DEF.\1j  LTS.I’.V.S. 
■Automatically  scaled  and  user-selected  clipped  examples  of  graphical  output  are  included  beneath 
the  code. 


PROGRAM  EASYPLT; 

USES 

GRAPH,  CRT. 

Coords,  Axes,  Mathlib,  Plotlib,  Form,  FloatDaf,  TypeDet,  Defaults; 
TYPE 


readfile  =  TEXT; 

ARRAY500  =  ARRAY [1 .. 500]  of  Float; 
VAR 

datafile 

Time,  Response 

i,  imin,  imax,  Ipts 

TSteps,  RSteps,  Tdecimals,  Rdecimals 

Tmin,  Tmax,  Rmin,  Raax:  Float; 

TWorld:  IITEGER; 


readf lie ; 
arraySOO ; 
IITEGER; 
IITEGER; 


ch  :  CHAR; 

BEQII  {Program} 

ASSIGKdataf  ile .  •  PLOT .  DAT  ■ ) ; 

RESETfdatafile) ; 
i  :=  0; 

WHILE  lOT  EOF(dataflle)  DO 
BEGII 

i  :=  i  +  1; 

READLKdatafile,  Timefi],  Responsefi]} ; 
EID; 

CLOSEfdataf ile) ; 

Ipts  i; 


Plotting  Routines  Folios 

} 

FindExtrema(Tiae,  Ipts,  Tmin,  Tmax,  imin,  imax); 
FindExtramaCRasponss,  Ipts,  Rmin,  Rmax.  imin,  imax); 
AxisScaleC 'TINE’ ,  TSteps,  Tmin,  Tmax,  TDecimals); 
AxisScalaC 'RESPOISE' ,  RSteps,  Rmin,  Rmax,  ROscimals); 
Def ineWorldCTVorld , tmin , Rmin , tmax , Rmax) ; 


Set  user-defined  defaults  and  the  size  of  the  plot  rectangle 

} 

SatDef suits; 

AxisColor  LightBlus; 

GEHStart ( ’ EASYPLT . GEM’ , Landscape ,0.68); 

ScrasnStart (0 . 8S ) ; 
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{ 

Draw,  label  and  naae  axes 

} 

LabellaxisClightgreen.TSteps,  Tdeciaals.  TVorld); 

Label YaxisCligbtgreen, RSteps .Rdeciula ,  TVorld) ; 
GemTextSettinga (Charter , 20) ; 
laneXaxisClightgraen,  'Tine  (seconds)*): 
8aBaYaxls(l'!ghtgrean,  'Response'); 

{ 

Plot  curve  data 

} 

PlotCarve(Yallov,  Ipts,  Linear,  Tine,  Response,  TVorld); 
ch  :=  ReadXey; 

PlotEnd; 

EHD. 


2.0  2.5  3.0  3.5  4.0  4.5  5.0  5.5  5.0  0.5  7.0  7.5  5.0 
lime  (secondt) 
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Appendix  II 

Test  Program:  MULTIPLT.PAS 


This  second  example  is  designed  to  illustrate  a  number  of  features  of  the  TP/Gem  interface.  Tiic 
code  is  a  significantly  modified  version  of  a  genuine  program  and  as  such,  the  actual  numbers  and 
plotted  curves  should  not  be  taken  too  seriously.  Examples  of  multiple  World  use,  global  variable 
changes,  user-defined  labelling,  NDC  dashed  lines  connecting  curves  from  different  Worlds  and 
logarithmic  axes  are  all  included  together  with  procedures  for  writing  screen  data  within  the  plot 
rectangle  and  an  application  of  the  unit  FORM  procedures. 


PROGRAM  NultiPlt; 

USES 

GRAPH. CRT. 

Coorda,  Axaa,  Mathlib,  PlotLib,  SBazier,  FloatDaf,  Fon,  TypaDat; 
TYPE 

ArraySOO  =  ARRAYfl . . 300]  OF  Float; 

VAR 
ch 

ynin,  yaaz,  xl,  yl,  z2,  y2,  Inin,  Inax,  Mnin,  NBax,  Tb, 

RqMia,  RqHax,  laabdal,  laabda2,  fq,  YqRaf,  Rqlora, 

RqMain,  RqMaax,  dlaabda,  IqFixad,  laabdaRat,  lanbdaFixad, 

RqRal,  RqFixad,  EflactivaCutof f ,  Cutoff Rat, 

CutoffFtxad,  Cl,  NCutoffFixad,  laabdaCutof f , 

HorkFunctlon 

i,  j,  Ipta,  IStapa,  MStapa,  lain,  iaax,  RqKStaps,  RqStaps, 

LogVozld,  RqWorld,  RqMVorld,  NVorld,  IDacinala,  NDaciaals, 
RqOaeiaala,  RqMOacinala 

lanbda,  laabdaT,  FracMlanbda,  Rq,  RqN,  «x,  ay 
HyLabala,  TlckPoan 
PROCEDURE  LotaOf Calculations; 

{ 

Tbs  cods  in  this  procsdurs  involvaa  calculation  of  array  data  that 
is  not  rslsvant  to  illustration  of  tha  operation  of  tbs  TP/GEN 
intarfaca.  Tha  coda  is  not  thsrafora  raplicatsd  aithin  this 
Appandix . 

> 

PROCEDURE  DisplayStuffOnScrsan; 

VAR 

TbStr,  LoasrStr,  UpparStr,  daltaStr,  EnargyStr,  ClStr, 

Rqstr,  CutoffStr,  fqStr,  lanbdaStr  :  STRUG; 

BEGII 

{ 

Title  plot  and  arits  data 

} 

SstColorfLightcyan) ; 

SstTaxtJnstifyfLaftTszt ,  BottoaTait); 

SstTaxtStyla(2,  HorizDir,  4); 

NoTtTo(DavX(ScraanDsTica,  0.0)  *  10,  DaaYfScraanDaaica,  1.0)  10); 

OntTsztf 'Background  T  ;  '}; 

Str(Tb;6:l,  TbStr); 

OutTazt(TbStr) ; 

OutTaxt('E'); 


:  CHAR; 


:  Float ; 


IRTEGER; 
Array 300; 
Array 100; 
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MoveTo(DaTX(ScrMiiO«vica,  0.0)  *  10.  DavYCScraanDavic*.  1. 0)^30); 
□utText( ’Ralaranca  Uavelangth  (1)  ;  '): 

Str(laBbdaK«l:5:3,  lubdaStr); 

□utTaztClaabdaStr) ; 

OutText ( ’  Bicrona ' ) ; 

NovaToCDavKScraaiiOavica,  0.0)  +  10,  DavYCScraanOavica,  1.0)  +  40); 
OutTaxtC’Ratarance  Raaponalvity  ;  ’): 

Str(RqRaf:8;6,  Rqatr); 

OutTaxt(Rqttr) ; 

MovaTo(DavI(ScraanOa*ica,  0.0)  *  10,  OavY(Scraai>Da*ica,  1.0)  +  50); 
OutTaxt( 'ElYactiva  Cutoll  Uavalangth  :  '): 

Str(CutoYfRal:S;3,  CutoflStr); 

□utTaxt(CutoffStr) ; 

OutTaxtC  ■icroBi’); 

NovaTo(DavX(ScraanDavica,  0.0)  *  10,  DavYCScraanOavica,  1.0)  >  60); 
OutTaxtClq  :  ■); 

StrCfqRaf ;8;6,  YqStr); 

□utTaxtCfqStr) ; 

NovaToCDavXCScraanOavica,  0.0)  *  10,  DavYCScraanOavica,  1.0)  *  80); 
OutTaxtC’RaXaranca  Vavalangth  C2)  ;  '); 

StrClaBbdaFixa<l:S;3,  laabdaStr) ; 

OutTaxtClaabdaStr) ; 

OutTaxt  C ‘  aicrona ' ) ; 

MovaToCDavXCScraanOavica,  0.0)  +  10,  DavYCScraanOavica,  1.0)  90); 

OutTaxt C’Ralaranca  Raaponaivity  :  '); 

StrCRqFixad:8:S,  Rqatr); 

OutTaztCRqatr) ; 

KovaToCOavXCScraanOaviea,  0.0)  +  10,  DavYCScraanOavica,  1.0)  +  100); 
OutTaztC’EfXactiva  Cutoll  Uavalangth  :  ’): 

StrCCutoffFixad:8;3,  CutoYfStr) ; 

OutTaztCCutollStr) ; 

OutTaxtC’  aicrona’): 

KovaToCOavXCScraanOaviea,  0.0)  +  10.  DavYCScraanOavica,  1.0)  ♦  110); 
OutTaxtC 'fq  :  ’); 

StrCfqFixad;8:6,  IqStr) ; 

OutTaztCfqStr) ; 

HovaToCDavXCSeraanOavica,  0.0)>  350,  DavYCScraanOavica,  1.0)  *  10); 
OutTaxtC ’Uork  Function  :  ’); 

Str CUorkFunction: 5 : 3 ,  EnargyStr) ; 

OutTaxtCEnargjStr) ; 

OutTaxtC’  a¥’): 

KovaToCOavXCScraanOaviea,  0.0)  *  350,  DavYCScraanOavica,  1.0)  *  20); 
OutTaxtC ’Cutoff  Uavalangth  :  ’); 

StrClaabdaCntoff :8:3,  Cutoff Str) ; 

OutTaztCCutoffStr) ; 

OutTaxtC ’  aicrona ’ ) ; 

KovaToCPavXCScraanOavica,  0.0)  *  350,  DavYCScraanOavica,  1.0)  *  30); 
OutTaxtC ’Cl  :  ’); 

StrCCl:S;3,  ClStr); 

OutTaxt CClStr); 
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PROCEDURE  SatValuas; 

BEGII 
Tb  :=  290; 

VorkFunction  0.219; 

Cl  ;=  0.267; 
laabdal  :=  1.0; 
laBbda2  :=  6.5; 
dlaabda  :=  0.02; 

EID; 

BEGII 

SatValnaa ; 

OpanFonC 'PtSi  Cutoff  avaluation’ ) : 

FloatItaH(Tb,  1,  ’Background  taaparatura  (K)  ;’); 

FloatItaa(VorkFuiiction,  3,  '  PtSi  Work  Function 

FloatltanCcl,  2,  *PtSi  raaponsivity  constant 

Floatltandanbdal,  1.  '  Louar  Spactral  Bound 

FloatItan(laabda2,  1,  '  Uppar  Spactral  Bound 

FloatItanCdlaabda,  2,  '  Spactral  Interval  :’); 

ClotaFom; 

FomControl  :=  1; 

WHILE  (FonControl  <>  100)  and  (FonControl  <>  101)  DO 
BEGII 

CASE  FonControl  OF 

1  :  GatRaalValue(Tb) ; 

2  ;  GatRaalValtta(HorkFunctlon) ; 

3  :  GatRaalValua(Cl); 

4  .-  GatRaalValue(laabdal) ; 

5  .-  GatRaalValua(lanbda2),' 

6  :  GatRaalValua(dlaBbda) ; 

EID;  {CASE} 

EID;  <WHILE> 

LotsOf Calculations ; 

FindEztraaadaabda,  Ipts,  Inin,  laaz,  ijiin,  inaz); 

FindEztranaCFracHlanbda,  Ipts,  Inin,  Hnaz,  inin,  inaz); 

FindEztraaa(Rq,  Ipts,  RqNin,  RqHaz,  iain,  laaz); 

FindEztraaa(RqN,  Ipts,  RqHain,  RqKaaz,  iain,  iaaz); 

AzisScals( 'HAVELEIGTI' ,  IStaps,  lain,  laaz,  IDsciaals); 

AzisScalaCFRACTIOIAL  EIITAICE',  MStaps,  Main,  Hnaz,  MDacinals); 

AzisScalaC ’PHOTOI  RESPOISIVITY’ ,  RqStaps,  RqMin,  RqNaz,  RqDacinals) ; 
yaaz  :=  RqHaz; 

LogScalaC 'LOG  TEST’,jain,  jraaz,  3);  {Only  included  as  an  ezajnple  of  use  of} 

{LogScala,  which  would  normally  be} 
{applied  in  association  with  a  log} 
{qualifiar  within  Procadura  PlotCurve} 

lain  2.6; 

laaz  6.0;  {Thasa  values  are  user-inspired,} 

Haaz  ;a  0.26;  {overwriting  the  auto-scaled  values} 

Rqlaz  ;•  6.0;  {calculated  within  AzisScale} 

RqHStaps  ;e  HStaps; 

RqHaaz  Ikaz; 

Daf inaHorld(RqWorld,  lain,  RqHln,  laaz,  RqHaz); 

DefiaeHorld(HVorld,  lain,  Hain,  Isuui,  Haaz); 

Daf ineWorld(RqHVorld,  lain,  RqHain,  laaz,  RqHaaz); 

Daf InaWorldCLogWorld,  lain,  yain,  laaz,  yaaz); 


WSRL-GD-02/91 


32 


{ 

S«t  tli«  *iz»  of  the  plot  roctangl* 

} 

GEMStart ( 'MULTIPLT . GEK ' . Landacapo ,0.68); 

Scr tenStart (0.85); 

Laballaxia  (■agauta,  IStapa.  IDacimala.  HqWorld); 

Side  ;=  InSide; 

AxiaTicka  :=  TwoSideaTicked; 

LogYAxia  (lightgreen.  LogWorld); 

Side  :=  OutSide; 

AxiaTicka  ;=  OneSideTieked; 

Ge«TextSettinga(Dutch,  20): 

laaaXaxia  (ligktmagenta,  'Wavelength  (Bicrona)’); 
laBaYaxia  (lightgreen,  ’LoglO(Photon  Reaponaivity) ’) : 

FOR  i:=  1  TO  S  DO 
BEGII 

MyLabelaCi)  :=  0.20  +  i*0.0S; 
yickPoanCil  •-  1 .240/KyLabelaCi3 ; 

EID; 

GaBTeitSettinga(S»iaa, 17); 

VarLabelTop  (lightgreen,  S,  TickPoan,  MyLabela,  2.  lain,  Inax); 
laaeTopAxiadightBagenta,  'Photon  Energy  (eV)'); 

PlotCurve  (blue,  Ipta,  Linear,  lambda,  FracMlaabda,  MWorld) : 
SBoothBe2ier(Yello«,  ipta,  0,  Linear,  lambda,  FracMlambda,  MWorld) 
SmoothBezierdightgreen,  Ipta,  0,  Linear,  lambda,  Rq, 
SmoothBezierdigbtred,  Ipta,  0,  Linear,  lambda,  RqH,  RqMWorld); 

{ 

Plot  daahed  linea 

> 

«DC(RqWorld,  CutolfFixed  ,  0.0,  xl,  yl); 

■DC(RqWorld,  CutoffFixed,  Rqlorm,  x2,  y2); 

PlotDaahedLine(Yello»,  Daahed,  Halfvidth,  xl,  yl,  x2,  y2) ; 
IDC(RqWorld,  Imin,  Rqlorm,  xl,  yl); 

PlotDaahedLine(Yello»,  Daahed,  Halfwidth,  xl,  yl,  x2,  y2); 
IDC(NWorld,  CutoffFixed,  MCutof fFixed,  xl,  yl): 

PlotDaahedLine (Yellow,  ShortDaahed,  Halfwidth,  xl,  yl,  *2.  y2) : 

D iaplayStuff Onscreen ; 
ch  :=  ReadXey; 

PlotEnd; 

Phown  Eiwnnr  (eV) 
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Appendix  III 

Test  Program:  GRAFTEST.PAS 


GRAFTEST  is  designed  to  demonstrate  the  application  of  the  procedure  SmoothBezier.  For  curves 
that  exhibit  high  curvatures  and  a  large  range  of  curvatures,  SmoothBezier  will  provide  a  sig¬ 
nificantly  better  fit  than  PlotBezier  at  the  expense  of  slower  operation  and  a  greater  memory 
requirement.  Note  that  the  choice  of  inside  ticks  with  unit  length  will  produce  a  complete  grid. 


{ 

Program  to  daaonatrata  plotting  routines.  Plots  frequency  response 
curves  of  the  fom: 

20I.OG  1.0  /  (1.0  -  (W/Wn)'2  +  2*ZETA«(W/Hn))  v  LDC(«/Wn) 
for  W/Wn  =  0.1. .10.0 

Separate  curves  for  ZETA  =  0. 1 ,0 .2.0.4,0.6, 1 .0 

} 

PROGRAM  GrafTast; 

USES  Graph,  Coords,  PlotLib,  SBezier,  Axes,  MathLib,  PloatOef,  TypeDef ; 

TYPE 

ZTabType  =  ARRAY [1..S]  OF  Float; 

COIST 

ZTab  :  ZTabType  =  (0. 1 ,0 . 2,0 .4,0.6, 1 .0) ; 

MPoints  =  80; 

VAR 

i,j,PlotVorld  :  IITEGER; 

Zeta,uu,dx,XMin,XKax  :  Float; 

X,Y,WorkX,«orkY  :  ARRAYCO . . 100]  OF  Float; 

U  :  ARRAY  Cl.. 2]  OF  Float;  {Complex  lumber} 

BEGII 

XHin  0.1; 

XMax  10. 0; 

Oaf ineVorldCPlotVorld , XMin , -40 , XMax , 20) ; 

{ 

Creates  a  new  world  referenced  by  “PlotVorld" 

World  is  bounded  by  (Xmin,  -40)  on  the  loser  left  comer  and 
by  (XMax,  20)  on  the  upper  right. 

) 

GemStart( ' GRAFTEST . GEM ' , Landscape, 0.7S) ; 

{ 

Enables  output  to  GEM  meta-file  GRAFTEST. GEN.  Orientation  is  leuidscape 
and  0.7S  of  the  plot  area  is  used  for  plotting  to. 

> 

GemWidth(B);  {Draws  GEM  lines  S  points  side) 

ScreenStart(0.75) ;  {Plots  to  0.7S  of  the  available  screen} 

TickLen  :=  1.0;  {length  of  ticks  on  tbs  axis} 

Side  Inside;  {plot  ticks  on  inside  of  plot  area  (normally  outside)} 
GaaTsztSattings (Charter,  12);  {Text  output  for  gem  is  Charter  font,  12  Point) 
LogXAxis(LightHagenta, PlotVorld);  {Draws  Z  Axis  in  light  magenta  colour  within 

World  defined  by  PlotVorld} 

LabelYAzis(LightHaganta, 6,0, PlotVorld) ;  {Draws  Y  Axis  in  light  magenta 

in  world  defined  by  PlotVorld. 

6  Ticks  and  numbers,  lo  decimals} 
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dx  (LoglOdMax)  -  LoglO(XMin))  /  IPoints; 

FOR  J  :=  1  TO  6  DO 
BEGIR 

Zeta  :=  ZTabCJ]; 

FOR  i  ;=  0  TO  RPointa  DO 
BEGII 

X[i]  :=  Pos«r(10,i  *  dx  *  LoglO(XMiD)) ; 

UClJ  :=  (1.0  -  X[iJ  • 

U[2]  :=  2.0  •  Zata  *  X[i]; 
uu  :=  U[l]  •  U[l]  +  UC2]  *  U[2]: 

U[l]  ;=  UCl]  /  uu; 

U[2]  :=  U[2]  /  uu; 

YCi]  :=  20  •  LoglO(SQRT(U[l]»U[l]  ♦  U[2]  •  U[2])); 
EHO; 

SmoothB«zier(Cyan,IIPoints>l ,  j.  LogLinear ,x,y,PlotWorld) ; 

{ 

Draas  a  amoothad  baziar  curva  in  PlotWorld 
j  sata  tha  lina  typa 
Log  X  axia,  Linaar  Y. 

X,  y  contain  data  pointa 

} 

PlotCurva (LightRad , MPointa . LogLinaar . x , y , PlotWorld) ; 

{ 

Draaa  tha  curva  through  x  and  y. 

Colour  ii  light  rad, 

MPointa  Bill  ba  plottad. 

Log  X  axia,  Linaar  Y  axia. 

} 

END; 

READLI : 

PlotEnd;  {Enda  acraan  graphica  and  cloaaa  GEH  mata-fila} 
WRITELK 'Monconvargance  ;  ’,MonConv): 

WRITELM(ErrorStr) ; 

END. 


■I  a  I 
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Appendix  IV 

Test  Program:  TESTDSM.PAS 


TESTDSM  provides  examples  of  the  use  of  the  discrete  smoothing  procedures  with  different  valiii's 
chosen  for  the  smoothing  parameter.  The  procedure  PlotMarkers  is  used  to  mark  the  individual 
points  defining  the  original  data  set. 


PROGRAM  DSMTust; 

USES  Graph, FloatDel,  TypaDat,  DSM. Coords. PlotLib.Axas.Bazler,  Crt; 
VAR 

*.y  .P.^^.d.ya  ,wgs  :  ARRAY  [0..101  OF  Float; 
u,v,dv,8x,sy  ;  ARRAYCO . . 100]  OF  Float; 
i.PlotWorld  :  IMTEGER; 
h,Rho  ;  Float; 

Ch  ;  chu; 

BEGIl 


-fAsaign  X  y  and  Wgs  arrays  initial  valuas} 


xCO] 

=  0.0;  y[0] 

*  3.0;  WgsfO] 

=  2.0 

xCl] 

=  1.0;  yCl] 

=  1.0;  WgsCl] 

=  1.0 

xC2] 

=  2.0;  y[2] 

=  2.0;  WgsC2] 

=  1.0 

X[3] 

’  3.0;  y[3] 

=  6.0;  WgsC3] 

=  1.5 

XC4] 

=  4.0;  yC4] 

=  5.5;  WgsC43 

=  1.0 

xCS] 

V  S.O;  y[S] 

=  6.0;  Wgs[S] 

=  1.0 

X[6] 

=  6.0;  yC6] 

=  7.0;  WgsCe] 

=  1.5 

XC7] 

=  7.0;  yC7] 

=  5.5;  WgsC7] 

=  1.5 

XC8] 

a  8.0;  yC8] 

»  4.0;  Wg8[8] 

=  1.0 

XC9] 

=  9.0;  yC9] 

=  3.5;  Wgs[9] 

*  1.0 

xClO]  :«  10.0;  yClO]  :=  1.0;  WgsflO]  :=  2.0; 

FOR  i  ;=  0  TO  100  DO 

u[i]  ; *  i  •  0.1;  -flnitialisa  u  array} 
h  :=  0.1; 

Rho  :=  0.00001; 

ScraanStart(0.7S);  {Start  scraan  graphics} 

GamStartf 'TESTDSM .GEM' .Landscape, 0 .65) ;{Start  Geo  graphics} 

DafinaWorldfPlotWorld, 0.0, 0.0, 10. 0,10.0) ;  {Define  «orld  PlotWorld) 

LabelXAxis (Green, 10,0, PlotVor Id) ;  {Assign  ticks  and  numbers  to  X  and  Y} 

LabalY Axis (Green, 10 , 0 , PlotUorld) ;  {axes} 

PlotMarkers(Blue, Star, 0.02, 11, Linear. x,y, PlotWorld);  {Plot  stars  as  markers} 

{at  data  points  for  x  and  y} 

laBeTopAxis(Maganta, 'FitDSM') ;  {Labels  top  axis  «ith  "FitDSM"} 

FOR  i  ;=  1  TO  6  DO  {five  Discrete  smoothing  curves  with  different  emphasis} 

{on  smoothing  (i=l)  to  curve  fitting  (i=5)} 

BEGIl 

Rho  ;=  Rho  a  lO;  {Rho  datermines  the  emphasis  -  see  manual  and  above} 
FitDSM(h,Rho,ll,z.y,Vgs,ys,b,c,d);  {compute  discrete  cubic  spline} 

FitCSdl ,  101  ,u,v ,x ,ys ,b,c.d)  ;  {Evaluate  cubic  spline  at  each  point} 

FitOSdl,  10i,u,dv,x,ys,b,c,d);{EvaluaCe  first  derivative  of  cubic} 

{spline  at  each  point  (not  plotted)} 

PlotCurva(i,  101,  Linear,  u,  v,  PlotWorld);  {Plot  the  nes  curve} 

EID; 
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PlotBezier(Yello9, 101 ,0 .Linear .u.v .wz.sy .PlotUorld) ; 

{ 

Plot  the  final  curve  using  beziera. 

Colour  is  Yellow. 

Plots  101  points.  Bo  dashes  (IDashes  =  0). 

Plot  is  against  a  linear  axis. 

Curve  is  plotted  through  points  defined  by  arrays  u  and  v. 
vx  and  vy  are  working  arrays.  These  are  used  by  PlotBezier. 
Curve  is  plotted  in  world  defined  by  PlotVorld. 

> 

While  not  KeyPreaaed  Do; 

Ch  ;=  RaadXey;  {clear  the  keyboard  buffer} 

PlotEnd; 

END. 


3T 
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Appendix  V 

Test  Program:  TESTMARKER.PAS 


PRQGRAH  T«stMark*r; 

USES  Graph. PlotLib, FloatD«l ,Typ«D«f.  Axes, Coords, Intarp; 

VAR 

x,y  :  Tlveetor; 

XX. yy  :  ARRAY[0. .200]  OF  Float; 
i.PlotWorld  :  IITEGER; 

Error  :  BYTE; 

BEGII 

FOR  i  :=  0  TO  10  DO 
BEGII 

xCi+1]  :=  i  •  0.1; 
y[i+l]  ;=  x[i+l]*xCi+l] ; 

EID; 

FOR  i  :=  0  TO  100  DO  xx[i]  ;=  i  •  0.01; 

CubicSplineFreedl  .x.y ,  lOl.xx.yy, Error) ; 

{ 

takes  11  data  points  Iron  arrays  x  and  y  and  returns  them 
as  101  data  points  in  arrays  xx  and  yy. 

> 

ScreenStart(0.7S) ; 

{ 

Initialise  screen  graphics  using  0.7S  of  the  screen  area 

} 

GenStart ( ' TESTHARK . GEM ’ . Portrait ,0.70); 

C 

Initialise  GEM  plotting.  Orientation  Portrait,  uses  0.7  of  the  plot 
area.  Gem  output  Bill  be  sent  to  file  TEST. GEM 

} 

Det ineWorld(PlotWorld,0 .0 ,0 . 0 , 1 .0, 1 .0) ; 

{ 

Define  a  world  indexed  by  PlotUorld. 
lower  left  is  (0,0),  upper  right  (1,  1) 

> 

LabelXAxisCMagenta, 10, 1 .PlotVorld) ; 

{ 

Draw  10  ticks  and  numbers  in  magenta  along  x  axis,  colour  is  magenta. 
Displays  1  decimal  point.  Drawn  in  world  defined  by  PlotVorld. 

} 

LabelYAxis (Magenta, 10, 1 .PlotVorld) ; 

{ 

Draw  10  ticks  and  numbers  in  magenta  along  y  axis,  colour  is  magenta. 
Displays  one  decimal  point.  Drawn  in  world  defined  by  PlotVorld. 

} 

PlotCurve(Cyeui, 11, Linear, xCl] .yCl] .PlotVorld); 

{ 

Draws  a  cyan  curve  of  x  vs  y  against  a  linear  axis.  Curve  is  drawn  in 
world  defined  by  PlotVorld. 

} 
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PlotHukersCLightHagcnta, Plus  ,0.02. 11  .Linear ,x[l]  ,y [1]  .PlotWorld) ; 

{ 

Draus  a  as  a  narkar  at  each  of  the  data  points  described  by  arrays 
X  and  y.  Markers  are  light  naganta,  drasn  in  sorld  PlotWorld  against 
a  linear  axis. 

Plots  11  Markers  at  scale  0.02  of  nomal  size  of  standard  device 

} 

PlotCurva(Green . 101 . Linear , xx , yy , PlotWorld) ; 

As  aith  above  PlotCurva  only  plots  xx  and  yy  arrays  which  sera  returned 
by  CubicSplinaFrae.  Careful  examination  shoas  the  difference  between 
the  smoothed  curve  and  the  plotted  original  data. 

} 

REAOLII ; 

PlotEnd;  {end  screen  and  gem  graphics} 

EHD. 
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Appendix  VI 

Test  Program:  SMOOTHTEST.PAS 


{»»»••«»»•*•*♦*»•**»*♦••*•***♦«•••*•*•*•*••••*•*•••*•••*•••••»»»*»*•*} 

{  PROGRAM  to  tost  the  snoothing  routine  Saoothlt.  } 


PROGRAM  SmoothTest; 

USES  Graph.Smooth.PlotLib,  Coords,  FloatDel.TypeDef ,  Axes; 

VAR 

x.y.z  ;  ARRAYCl. .1024]  OF  Float; 
dx  ;  Float ; 

i.PlotVorld  :  IITEGER; 

BEGII 

RAIDOMIZE;  {Seed  randon  number  generator} 
dx  ;=  0.01; 

FOR  i  :=  0  TO  1000  DO 
BEGII 

xCitl]  :=  i  *  dx; 

z[i+l]  :=  SII(x[i+l]);  {sine  values,  x  versus  z} 

yCi+1]  :=  zCitl]  *  0.5  •  (Random-0.5);  {sine  values  from  z  with  random} 

{noise  added} 


£10; 

GEMStart ( ■ SMOTEST . GEM ' . Landscape ,0.68); 

ScreenStart(0.76) ;  {Start  Screen  Graphics} 
DofineWorld(PlotWorld,0,-2,10,2);  {Define  a  world  to  plot  in} 
LabelXAxisCRed, 10,0,PlotVorld) ;  {Put  numbers  and  ticks  on  X  axis} 

LabelYAxis(Red,4,0,PlotMorld) ;  {Put  numbers  emd  ticks  on  Y  axis} 

PlotCurve(M2:genta,  1001, Linear  ,x,z,PlotUorld) ;  {Plot  sine  curve  x  vs  z} 
PlotCurveCRed, 1001, Linear, x.y.PlotWorld);  {Plot  sine  curve  with  "noise"} 
SmoothftCy, 1001, 100); 

{ 

takes  array  y  for  smoothing, 
smooths  1001  points, 
smooths  over  100  points. 

} 

PlotCurve(Green, 1001 , Linear ,x,y ,PlotWorld) ;  {Plots  the  smoothed  curve} 

{returned  in  y} 

READLI; 

PlotEnd;  {End  screen  graphics} 

EID. 
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Appendix  VII 

Test  Program:  FORMTEST.PAS 


PROGRAM  FonnTest; 

USES 

Crt ,  Fora ,  FloatDal ; 

VAR 

SomeVunber  ;  IITEGER ; 

Lanbdal,  Laabda2,  DLaabda, 

Tb,  CutOn,  UorkFuaction,  Cl  :  Float; 
astring  :  String; 

PROCEDURE  Sat Values; 

BEGII 
Tb 

WorkFunction 
Cl 

CutOn 
Lambda 1 
Lambda2 
DLambda 
Somalumber 
astring 
EKD; 

PROCEDURE  DisplayForm; 

{ 

This  Procedure  sill  sat  up  the  screen  display  for  a  desired  input  fora. 
Procedures  used  axe  - 

OpenFora  *  Will  display  at  the  top  of  the  screen  the  desired  heading. 

The  heading  is  passed  as  a  parameter  within  quotes. 
AdvancsLine  -  Allows  a  blank  line  to  be  inserted  on  screen. 

CloseFom  -  Displays  a  footer  for  the  input  fora.  It  displays  the 
message  ’Press  PgDn  OR  ESC  to  Accept  Values'. 

Float Item, 

Integerltea,  Used  for  placing  prompts  on  screen,  Stating  what  vuiables 
Stringltem  are  associated  with  each  prompt  and  the  number  of  decimal 
points  to  be  displayed  with  each  veuriable 

-  Ceill  procedure  Floatitem,  Integerltem  or  Stringltem  depending 
upon  what  sort  of  data  item  is  on  a  given  line. 

-  Parameters  to  the  function  calls  are  (in  order) 

1.  Variable  to  store  the  input  data. 

2.  lumber  of  Decimal  Points  to  display.  For  Integer  and  String 
Values  this  number  is  ignored. 

3.  The  prompt  to  appear  for  the  value  at  a  given  line. 

Mote  that  the  prompt  will  commence  in  Column  5.  The  cursor  will 
be  placed  one  space  after  the  end  of  the  prompt . 

A  form  should  be  built  up  in  the  following  manner 

1  -  OpenFormC ’Form  Heading’); 

2  -  Floatitem 
or  Integerltem 
or  Stringltem 


2S0; 

0.213429; 

-0.267E-09: 

3.4; 

1.0; 

6.S; 

0.02; 

8; 

’Hello  World’; 
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Thas*  should  appear  on  each  line  uhars  thara  sill  bs 
an  input  proapt.  Thay  aay  ba  intarsparsad  aith  a  call  to 
tha  AdvancaLina  procadura  in  order  to  iaprova  readability 
or  group  ralatad  itau  (or  both) .  Thay  should  ba  in  tha 
order  data  aould  noraally  ba  antarad  at  tha  pronpts 

3  '  ClosaFon: 

4  -  Having  built  up  tha  screen  it  ia  nacaasary  to  add  tha  coda 

to  read  tha  data  valuas. 

To  do  this  coda  aust  ba  insartad  batsaan  tha  beginning  and  end 

of  tha  Casa  Statanant. 

It  takas  tha  following  format  ; 

ForaControl  :=  1;  or  shichavar  field  cursor  goes  in  first 

Casa  ForaControl  of 

-  Start  inserting  code  hare  as  naedad  - 

1  ;  CatRaalValua(Tb); 

2  :  GatIntagerValue(C2) ; 


■  ;  GatStringValua(InputString) ; 

-  after  the  last  prompt  no  further  code  is  required  - 

End; 

Each  proapt  aust  ba  assigned  a  number  vhich  corresponds  to 
tha  position  of  tha  proapt  in  the  fora  layout  defined  above. 

First  Proapt  is  1 
Second  is  2 

and  so  on.  Tha  nuaber  MUST  ba  follouad  by  tha  character, 
note  that  tha  ‘Gat’  procaduras  accept  parameters  vhich 
corrtapond  to  tha  variables  defined  above 
S  -  The  fora  should  nos  ba  coaplata.  To  run  it,  call  DisplayFora  from 
your  application. 


BEGII 

OpanForaC'  PtSi  CUTOFF  EVALUATIOI’) 
FloatltaaCTb,  1,  * 

FloatltaaCVorkFunction,  3,  ' 

FloatltaaCCl ,  2,  ' 

AdvancaLina; 

FloatltaaCCutOn,  1,  ' 

FloatlteaCLaabdal ,  1 ,  ’ 

FloatItan(LaBbda2,  1,  ’ 

FloatltaaCDLaabda,  2,  ’ 

AdvancaLina; 

IntegerlteaCSoneluBber ,  2,  ' 

AdvancaLina; 

StringItaaCastring,  0,  ’ 

ClossFora; 


Background  Tenparature  (K) :  ' ) ; 
PtSi  Work  Function  (aV) 

PtSi  Responsivity  Constant  Cl 

Cuton  Wavelength  (microns)  . ' ) ; 
Loser  Spectral  Bound  (microns)  :'); 
Upper  Spectral  Bound  (aicrons)  :'); 
Required  Step  Size  (microns)  :'); 

Some  lumber  (<  10)  :'); 

Enter  a  String  hare  : ' ) ; 
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For»Control  :=  1; 

While  (ForaControl  <>  100)  and  (FornControl  <>  101)  Do 
Begin 

Case  FoxnControl  of 

1  :  GetRealValue(tb) ; 

2  ;  GetRealValue(WorkFunction) : 

3  :  GetRealValue(Cl) : 

4  ;  GetRealValue(CutOn) : 

E  :  GetRealValue(Laabdal) ; 

6  ;  GetRealValtte(Laabda2) : 

7  ;  GetRealValne(DLaabda) ; 

8  ;  Begin  {Eianple  of  validity  checking} 

GetIntegerValue(SoBeIu>ber) ; 

IF  Soaelnaber  >=  10  TBEI  {Value  out  of  range} 

BEGII 

GoToXY(10.24); 

WriteC'This  number  should  be  less  than  10*); 

ForaControl  :=  8; 

EID 

ELSE  {Ensure  Error  message  disappears} 

BEGII 

GoToXY(10.24) ; 

WriteC' 

EID; 

End; 

9  :  GetStringValue(astring) ; 

EID;  {Case} 

End; 

EID; 

BEGII 

SetValues; 

DisplayForm; 
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